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Abstract 

It has recently been suggested that the presence of a plenitude of light axions, an Ax- 
iverse, is evidence for the extra dimensions of string theory. Wc discuss the observational 
consequences of these axions on astrophysical black holes through the Penrose superradiancc 
process. When an axion Compton wavelength is comparable to the size of a black hole, the 
axion binds to the black hole "nucleus" forming a gravitational atom in the sky. The occupa- 
tion number of superradiant atomic levels, fed by the energy and angular momentum of the 
black hole, grows exponentially. The black hole spins down and an axion Bose-Einstein con- 
densate cloud forms around it. When the attractive axion self-interactions become stronger 
than the gravitational binding energy, the axion cloud collapses, a phenomenon known in 
condensed matter physics as "Bosenova" . The existence of axions is first diagnosed by gaps 
in the mass vs spin plot of astrophysical black holes. For young black holes the allowed values 
of spin are quantized, giving rise to "Regge trajectories" inside the gap region. The axion 
cloud can also be observed directly either through precision mapping of the near horizon 
geometry or through gravitational waves coming from the Bosenova explosion, as well as 
axion transitions and annihilations in the gravitational atom. Our estimates suggest that 
these signals are detectable in upcoming experiments, such as Advanced LIGO, AGIS, and 
LISA. Current black hole spin measurements imply an upper bound on the QCD axion de- 
cay constant of 2 ■ 10^^ GeV, while Advanced LIGO can detect signals from a QCD axion 
cloud with a decay constant as low as the GUT scale. We finally discuss the possibility of 
observing the 7-rays associated with the Bosenova explosion and, perhaps, the radio waves 
from axion-to-photon conversion for the QCD axion. 
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1 Introduction and summary 

Black holes are among the most fascinating systems in astrophysics, and the most mysterious 
objects in quantum gravity and string theory, for a long time serving as principal "theoretical 
laboratories" for exploring non-perturbative gravitational dynamics. The purpose of this paper is 
to initiate a detailed study of the exciting possibility that astrophysical black holes may serve 
as actual laboratories for the discovery of new elementary particles. 

There are several reasons why we believe this possibility is realistic. On a purely phenomeno- 
logical side, black hole observations are routine practice in nowaday astronomy (see, e.g., [2] for 
a review). About 40 stellar mass black holes in X-ray binaries in the Milky Way and neighboring 
galaxies have been identified with masses in the range ~ 5 ^ 20Mq. Supermassive black holes 
with masses ~ 10^ -i- 10^^ Mq have been found in centers of many galaxies including the Milky 
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Way and believed to be hosted by nearly all of the galaxies. Also, the first intermediate mass 
(~ 100 -i- 10^ ^0) candidates have been identified. 

Following the evolution of binary systems or measuring the velocity dispersion of stars rotating 
around galactic centers allows to determine black hole masses. Most crucially for what follows, 
recent advances in X-ray astronomy and in numerical magnetohydrodynamical simulations of the 
accreting gas in the Kerr metric open the possibility for a detailed exploration of the near-horizon 
region and, as a consequence, for high precision neasurements of black hole spins [3l H]. First 
estimates for the angular momentum of several black holes have already been delivered [S], often 
suggesting high values for the spin, although at the moment different techniques sometimes give 
rise to conflicting results [6] . 

In the future, apart from improvements of traditional astronomical techniques for observing 
the near horizon environment and its better theoretical modeling, a unique probe of the black 
hole geometry will be provided by low frequency gravitational waves observatories, such as LISA 
[7] or AGIS, a gravitational wave detector based on atom interferometric techniques P H]. For 
the purpose of testing the near horizon geometry the most promising candidates are the so-called 
extreme mass ratio inspirals — stellar mass compact objects captured by supermassive black hole 
in the galactic center (see, e.g., [10]). LISA and AGIS are expected to detect about a hundred of 
such events per year. Each such measurement allows not only to determine the black hole spin 
and mass with an exquisite accuracy, 10~^ -ir 10~^, depending on the details of a particular event, 
but also to check whether higher order metric moments, up to 6 4- 7, agree with their values for 
the Kerr geometry. 

This ongoing observational progress indicates that we are witnessing the dawn of precision 
black hole physics. Undoubtedly, black hole observations will be of great value for astrophysics, 
however it is natural to inquire whether these data may be useful for beyond the Standard Model 
physics as well, given that it will provide a rare test of non-linear gravity. However, possibly 
contrary to naive expectation, it turns out quite challenging to find modified gravity theories 
which would predict deviations from general relativity near astrophysical black holes and would 
not contradict current gravity tests. One candidate class of modified theories of gravity affecting 
black hole dynamics are models of Higgs phases of gravity, where black hole no-hair theorems can 
be violated [TT]. 

In this paper we explore a less exotic possibility to test fundamental physics with precision 
black hole observations. It is related to the famous Penrose process, a mechanism to extract 
energy and angular momentum from rotating black holes [121 Il3|. As reviewed in detail below, 
this process, known as superradiance, when applied to waves rather than particles [HI dSl ITB] . 
gives rise to a spin-down instability of a rotating black hole [TTJ dSl [T21 [20], if a massive boson with 
a Compton wavelength of order the black hole gravitational radius is present in nature. As we 
will see, this instability turns rotating astrophysical black holes into sensitive detectors of bosons 
with masses in the range fi ~ 10^^ -7- 10^^^ eV. Before focusing on the observational consequences 
of the superradiant instability, let us review why it is natural to expect ultra-light bosons in the 
theory that transform astrophysical black holes in probes of fundamental physics. 

A natural situation giving rise to a particle of a small, but non-vanishing mass is when this 
particle is a (pseudo)Goldstone boson of a spontaneously broken global symmetry, which is also 
explicitly broken by non-perturbative effects. Probably the best motivated candidate for such 
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a particle is the QCD axion (f)a — a pseudoscalar particle coupled to the QCD instanton number 
density via 



Note that at the classical level Se is invariant under the Peccei-Quinn (PQ) symmetry (pa 
(pa + const, so that the QCD axion is indeed a (pseudo)Goldstone boson with fa being the scale 
of spontaneous symmetry breaking. This symmetry is explicitly broken by the QCD instanton 
effects that generate the axion potential giving rise to a solution for the strong CP problem — the 
primary motivation for the QCD axion. As a result the QCD axion acquires a mass equal to 



The Compton wavelength of the QCD axion with a high symmetry breaking scale ^ 10^^ GeV 
matches the size of stellar mass black holes and, consequently, can affect their dynamics, suggesting 
that this part of the parameter space for the QCD axion can be explored through black hole 
observations. 

There are several reasons why this conclusion is very important. First, non-gravitational 
interactions of the QCD axion with the rest of the Standard Model particles are very suppressed at 
these high values of fa- As a result this part of the parameter space can not be easily probed by any 
other means, either laboratory or astrophysical. Second, in many "generic" string constructions, 
i.e., in compactifications where the extra- dimensional manifold is neither highly anisotropic, nor 
strongly warped, the values of fa are naturally around the grand unification scale Mqut — 2 x 
10^^ GeV [21]. Finally, as elaborated in more detail in section [sj finding the QCD axion with 
fa ~ McuT would indicate that the baryon-to-dark matter ratio varies on length scales longer 
than the observed part of the Universe and its local value is determined by anthropic reasoning. 
Discovery of the QCD axion in this regime would be further evidence for enviromental selection 
already suggested by the cosmological constant problem, and by the string landscape. 

There is an even stronger and more direct link between the QCD axion and the landscape of 
string vacua, a link that gives rise to the expectation of a plenitude of light axion-like particles, an 
axiverse pj| - this same link also suggests the existence of many massless vectors, whose massive 
superpartners may be discovered at the LHC |22] • In string constructions, an axion usually arises 
as a Kaluza-Klein (KK) zero mode of a higher-dimensional antisymmetric form field. Such zero 
modes have a purely topological origin: they are labeled by non-contractable cycles in the extra- 
dimensional manifold. Non-contractable cycles allow for non-trivial gauge field configurations with 
a vanishing field strength, the so called Wilson lines. These configurations do not carry energy 
and correspond to zero KK modes at the perturbative level. They only acquire a mass due to 
non-perturbative effects. 

Interestingly, the very same ingredients that give rise to the string axiverse, higher-dimensional 
form fields and non-trivial cycles in the compactification manifold allowing also to turn on gauge 
fluxes, also give rise to the string landscape of 10^°*^ or so vacua. In order to allow for the tuning of 
the cosmological constant at the ~ 10^^^° level, as required by observations, the compactification 
manifold should contain of order few hundred cycles, given that the total amount of flux quanta 
for a cycle is typically limited by a number around ten in order to stay in a perturbative regime. 




(1) 




(2) 
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Figure 1: Axionic Black Hole Atom: The spinning black hole "feeds" superradiant states form- 
ing an axion Bose-Einstein condensate. The resulting bosonic atom will emit gravitons through 
axion transitions between levels and annihilations and will emit axions as a consequence of self- 
interactions in the axion field. 



Consequently, one may expect hundreds of axion-like particles in a given string compactification. 
However, a plenitude of cycles does not yet guarantee the presence of a plenitude of axions. There 
is a number of effects in string theory that could produce a large axion mass, such as branes 
wrapping the cycles, and fluxes. One can roughly estimate the number of light axions as being 
determined by the number of cycles without fluxes — presumably, around one tenth of the total 
number of cycles. Still this leaves us with the expectation of several tens of axion-like particles. 

The discovery of a plenitude of particles in our vacuum with similar properties but different 
masses supports the idea of a plenitude of vacua, as both the axiverse and the multiverse are 
dynamical consequences of the same fundamental ingredients. 

The masses of string axions are exponentially sensitive to the sizes of the corresponding cycles, 
so one expects them to be homogeneously distributed on the logarithmic scale. However, given 
that the QCD ^-parameter is constrained to be less than 10"^*^, non-perturbative string corrections 
to the QCD axion potential should be at least ten orders of magnitude suppressed as compared 
to the QCD generated potential. It is then natural to expect many of the axions to be much 
lighter than the QCD axion; these are the axions whose mass is dominated only by these small 
non-perturbative string effects. 

The implicit, and very plausible assumption behind this line of reasoning is that there is no 
anthropic reason for the existence and properties of the QCD axion. Consequently, these properties 
should follow from the dynamics of the compactification manifold, rather than being a result of 
fine-tuning, and the QCD axion should be a typical representative among other axion-like fields. 
A priori we expect tens (or even hundreds) of light axions, it would be really surprising if the 
QCD axion turned out to be the single one. 
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These arguments motivate us to look not only for the QCD axion, but for axions in the entire 
mass range /i ~ 10~^ -7- 10~^^ eV, where they can affect stellar or galactic astrophysical black 
holes through superradiance. Let us summarize now the major features of superradiance and its 
principal observational consequences. 

Superradiance [m [151 IB] is the phenomenon of wave amplification during scattering off a 
rotating black hole which takes place whenever the wave frequency u and the magnetic quantum 
number m satisfy the superradiance condition 

< u < mw+ , (3) 

where w+ is the angular velocity of the black hole horizon defined as 

w^ = - (4) 

2r, 1 + Vl - {a/r,y 

with Tg being the gravitational radius of the black hole and < a < is the spin-to-mass ratio. 
Superradiant amplification may lead to an instability under certain circumstances. One example, 
admittedly not a very practical one, is the "black hole bomb" by Press and Teukolsky [2^ a 
rotating black hole surrounded by a spherical mirror. A single photon introduced in the system, or 
created by quantum fluctuations, with quantum numbers satisfying the superradiance condition 
^ gives rise to an exponentially growing number of photons inside the mirror through a chain of 
consequent amplifications at the black hole horizon and reflections from the mirror. 

Remarkably, nature provides such a mirror in the presence of a massive boson [T71 [18], [191 120] . 
Massive bosons, in our case axions, have bound Keplerian levels in the gravitational field of a black 
hole. This allows for a black hole to release its spin by populating levels satisfying the condition 
([3]), see Fig. [l} This creates an axionic Bose-Einstein condensate (BEG) cloud rotating around 
the black hole. The process is only efficient if the Compton wavelength of an axion is comparable 
to the black hole size. As a result, the energy spectrum of superradiant levels is quantized and 
very close to the spectrum of a hydrogen atom. The superradiant axion cloud loses its energy 
and momentum by gravitational wave emission associated with axion transitions between different 
"atomic" levels, and with axion annihilations to gravitons. Another important loss mechanism is 
related to non-linearities in the axion potential and results in the emission of axions. Finally, the 
whole system may be fueled by energy and spin inflow from the matter accreting onto the black 
hole. 

In Fig. |2]we show the region of black hole mass and axion mass parameter space that is affected 
by superradiance. The dark area outlines the region where a superradiant cloud has enough time 
to be built up during the lifetime of the Universe for a maximally spinning black hole. The lighter 
region shows the part of the parameter space where the superradiant spindown rate is faster than 
the spinup rate due to Eddington accretion, so that Eddington accreting black holes in this mass 
range lose their spin as the cloud develops. 

There are three major classes of observational signatures associated with the axion cloud that 
we are going to discuss in the current paper: 

• Gaps in the black hole "Regge plot" (spin vs mass plane) — the absence of rapidly rotating 
black holes when their gravitation radius matches the Gompton wavelength of an axion. 
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Figure 2: The part of the black hole and axion parameter space potentially affected by super- 
radiance. For axion and black hole masses in the colored region the time required to create a 
substantial axion cloud is shorter than the age of the Universe. For masses in a light colored 
region the superradiance rate is faster than the Eddington accretion rate. 



• Direct gravitational wave signals from the rotating axion cloud. For the QCD axion this 
signal falls into the sensitivity band of the Advanced LIGO interferometer. 

• Modification of the near-horizon metric due to the presence of the axion cloud. 

The main goal of this paper is to evaluate the observability of these classes of signatures in 
near future experiments. Detailed quantitative predictions for gravitational wave emission rates 
and waveforms, as well as the change of the templates for the near horizon metric is likely to 
require numerical work. The reason is that the axion BEG cloud is a rich and complicated 
dynamical system with many processes that need to be taken into account in order to have an 
accurate description of its behavior at cosmological time scales. One important process that poses 
a challenge for an analytical treatment is the effect of axion self-interactions on superradiance. 
Here we limit ourself to a qualitative, semi-analytical discussion of the major physical processes 
involved. Our analysis indicates that depending on the relation between the axion mass and the 
black hole size all three classes of signatures can be observable for different systems. This provides 
strong motivation for further numerical analysis of axionic superradiance. 

Let us briefly summarize our main results. In Fig. |3]we show the black hole Regge plot for two 
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Figure 3: The regions in the black hole Regge plane affected by superradiance for the QCD axion 
with m = 3 ■ 10~^^ eV (the upper panel) and for a lighter axion with m = 10~^^ eV. The data 
points correspond to spin measurements obtained by fitting the thermal continuum X-ray spectra 
[5] . Old black holes are expected to be found in the shaded region, where they are not affected by 
superradiance. Young black holes may be found also on the dashed colored lines inside the gap. 
Different colors correspond to superradiant levels with different orbital quantum numbers /. 
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different axion masses. The upper plot corresponds to an axion mass of /Xa = 3 ■ 10^^^ eV, which 
is the mass of the QCD axion for = 2 ■ 10^'' GeV, and the lower illustrates the effect of a lighter 
axion with /Iq = 10~^^ eV. Black holes in the shadowed region are not affected by superradiance 
during the age of the Universe. Sufficiently old black holes are expected to be found in this region. 
In the plot corresponding to the QCD axion we also present existing spin measurements; this data 
suggest an upper bound on the QCD axion decay constant of 2-10^^ GeV. 

Note that the uncolored gap regions exhibit internal structure (dashed lines), reflecting the 
quantized behavior of the gravitational black hole atom. Black holes may stay on these lines, 
the "Regge trajectories", for cosmological time scales. In order to understand this behavior, we 
should note that there is a number of different superradiant bound states labeled by different 
angular quantum number /. The instability time-scale rapidly drops down as the orbital quantum 
number / increases. Therefore, the spin-down process of the black hole is initially driven by the 
level with the minimum value of / for which the superradiance condition ^ is satisfied, and stops 
when enough spin has been extracted, so that the superradiance instability rate goes to zero — 
the superradiance condition is now saturated. At this point the black hole finds itself on one of 
the dashed lines of the Regge plot. 

Next, we should expect that the black hole spin-down proceeds with the growth of the level 
with the next-to-largest superradiance rate, i.e., the / -|- 1 level. However, further spin-down is 
damped when the black hole reaches one of the Regge trajectories as a consequence of axion self- 
interactions. The axions bound to the level that is no longer superradiating serve as an axially 
asymmetric perturbation to the system that mixes superradiant with non-sup err adiant levels and 
shuts off the black hole spin-down process. This is similar to introducing a non-spherical defect 
on the mirror in the Press- Teukolsky black hole bomb: photons reflected from the mirror exit the 
superradiant region Q and are now absorbed by the black hole, turning the instability off. As a 
result of these axion non-linearities, black holes stay on the Regge trajectories until various axion 
loss processes, or a violent accretion event, dissipate the cloud down to a small enough size, such 
that non-linearities no longer inhibit the instability. Then the black hole rapidly jumps to the 
next Regge trajectory by populating the level with a larger orbital quantum number. 

Axion non-linearities trigger yet another dramatic effect during transitions between Regge 
trajectories. Namely, attractive axion self-interactions result in a catastrophic instability of the 
axion cloud and its subsequent collapse. The analogous effect has been observed in laboratory 
BEC systems with attractive self- interactions and is known under the name "Bosenova" [25] . We 
find that, depending on the parameters, transitions between Regge trajectories proceed through 
a series of tens to hundreds of Bosenova events. These events produce gravitational waves, and, 
in the case of the QCD axion, gamma- and X-rays that may be detectable on the Earth. 

The black hole may also emit gravitational waves at the observable level due to axion transitions 
between different levels and annihilations to gravitons. Accurate prediction of the strength and 
duration of the signal requires a more detailed analysis of the dynamics; in the current paper we 
limit ourselfs to qualitative estimates which imply that this signal is detectable. It is especially 
exciting that the gravitational wave signal from the QCD axion cloud around stellar mass black 
holes falls into the Advanced LIGO frequency band and turns it into a particle discovery machine. 
In Fig. |4]we present the estimated signal strength as a function of the axion and black hole masses 
for a source at a 20 Mpc distance from the Earth. We choose a coherent integration time for the 
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Figure 4: The contour plot of constant gravitational wave signal from axion transitions between 
the Qg and the 5g levels for a black hole located at 20 Mpc away from the Earth. The projected 
sensitivity curves of Advanced LIGO [2S] and Einstein Telescope [2Z] assume 10^ seconds of a 
coherent integration time. 



signal of 10'^ second , when drawing the sensitivity curves. 

The paper is organized as follows. We start with explaining the basics of superradiance in 
Section [2j We argue that the non-relativistic approximation is accurate for describing superradiant 
levels for most of the parameter space, and adequate at least qualitatively at all values of the 
parameters. We present semi-analytical results for the superradiant rates, that are accurate over 
a large part of the parameter space and enough for our purposes. 

In Section |3] we discuss the processes that determine the dynamics of the superradiant cloud — 
gravitational wave emission and axion non-linearities. We provide approximate expressions for the 
relevant transition rates, and for characteristic timescales and masses of the superradiant cloud 
when different processes, such as Bosenova collapse, happen. 

In Section|4]we combine all of the above ingredients and discuss how the superradiant instability 
develops and what are the associated observational signatures. We discuss in detail the black hole 
spin-down and of the Regge trajectories. We then estimate the gravitational wave signals during 
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transitions between Regge trajectories and briefly comment on prospects for directly detecting tlie 
axionic cloud around supermassive black holes during extreme mass ratio inspirals. 

We also focus on the potential reach for the QCD axion, and some possible QCD axion specific 
signatures related to its direct coupling to Standard Model fields. A particularly intriguing effect of 
the QCD axion is that the effective value of the QCD ^-parameter may become of order one inside 
the cloud. This possibly destabilizes nuclei in the accretion disk and results in 7-ray signals and 
exotic X-ray lines from the nuclear decay products in the black hole vicinity. We finally entertain 
the possibility of radio waves from axion-to-photon conversion in the near-horizon magnetic field. 

Section [5] elaborates on an issue which is aside from the main line of the paper, but still a very 
important part of the theoretical motivations for string axions in the mass range probed by the 
black hole superradiance. Namely, all these axions are "anthropic" — their initial misalignment 
angle needs to be tuned to an atypically small value in the observed part of the Universe. It has 
been known for a long time [2H], that this is not a problem for a single QCD axion, and here we 
discuss what changes if several anthropic axions are present. We conclude in Section [6] 



2 Spectroscopy of superradiance 



In this section we review the spectroscopy of superradiant levels around a rotating black hole. 
Throughout the paper we are using the Boyer-Lindquist coordinates for the spinning black hole 
metric |29l 



ds^ = -(1 ^)dt^ ^ dtdd) + —dr^ + T.d9^ + ^ sin^ 

J 



S = -|- cos^ 6^ , A = (r — r+)(r — r_) r± = Vg ± \ — , = TT? t'g = G^M , (5) 

V M 

where M and J are black hole's mass and spin, respectively. The physical horizon corresponds to 
the larger root of A, which is r = r+. 

One of the most fundamental equations for the understanding of superradiance is condition 
(|3]), so let us review how it arises, following [29]. Interestingly, one has to know remarkably little 
about the Kerr geometry to derive ([s]). The first fact one needs is that the black hole metric in the 
Boyer-Lindquist coordinates is invariant under time translations, or more formally it possesses 
the Killing vector 'H'^d^ = dt- Second, it possesses another Killing vector related to rotational 
invariance, J^d^ = d^p and the linear combination 

g = n + w+j, 

is normal to the horizon and null there. In the above equation is given by Q. 
Now, let us consider an incoming wave of a scalar field of the form 

(j) = e-''^'+'"'^f{r,9) + h.c. . 

The conserved energy flux of this field is given by 
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where and C are the energy-momentum tensor and the Lagrangian density for 0, respectively. 
For a space-time region between two constant time shces, conservation of the current imphes 
that the time-averaged energy flux at the infinity is equal to the time-averaged energy flux through 
the black hole horizon. The latter is equal to 

{P^G^ = -{{dt(p + io+d^(p)dt(P)=uj{io-mio+)\f\^ , (6) 

where the ^f^i^-term in the energy-momentum tensor dropped out because vectors Q and H are 
perpendicular at the horizon. We see that the energy flux is negative in the superradiant frequency 
range indicating that the wave gets amplified in this regime. The argument changes a bit when 
the frequency u corresponds to the discrete spectrum so that the energy flux at the infinity 
necessarily vanishes. In this case the only way to reconcile the flux (j6]) at the horizon with the 
energy conservation is for the frequency u to acquire an imaginary part, so that the time- averaged 
energies on the two constant time slices are not equal any longer. For the real part of the frequency 
in the superradiant interval ([s]) the imaginary part should be positive indicating the presence of 
an instability. 

One important consequence of the instability condition ^ is that the superradiant levels are 
always in a (quasi)non-relativistic Keplerian regime. Indeed, in this regime the real part of the 
frequency follows the hydrogen spectrum 

, (7) 

where n = n + Z + 1 is the principal quantum number, / is the orbital moment and a = ^a'^g. For 
such a level the velocity of the particle is 

- . (8) 
n 

On the other hand, if we approximate the frequency in the superradiance condition ([3]) by the 
axion mass, a;^ ~ /^a? the condition translates into a bound 

a<mw+ = ——, (9) 

where we made use of the expression Q for We see now that the velocity for superradiant 
states may be at most moderately relativistic, 

^ 1 m a 1 , , 

2 n r_|_ 2 

where the bound is saturated at the upper boundary of the superradiant range ^ for extremal 
black holes a = r+, at n = and large / = m S> 1. 

In principle, this argument does not exclude the presence of a family of non-hydrogenic unstable 
bound states, however numerical results of |30] confirm that all superradiant states are hydrogenic. 

We can also estimate the size of the axion cloud as 

Tc ~ , (11) 
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which is always significantly larger than the black hole gravitational radius, as a consequence of 
(§■ 

These estimates provide the following physical picture of the superradiant axion cloud. The 
cloud is composed of a wave packet of the axion field rotating on a Keplerian orbit around the 
black hole. This axion wave packet always has a tail that goes into the near-horizon ergosphere 
region and gets amplified there leading to the exponential growth of the number of axions in the 
packet. 

Given the complexity of the Kerr metric it is not surprising that a precise analytical expression 
for superradiant rates is unavailable (partial numerical results can be found in [SO])- However, 
the above physical picture gives rise to two useful analytical approximations for the superradiant 
rates. Before introducing them let us recall that the massive Klein-Gordon equation in the Kerr 
background allows separation of variables [31] with the following simple ansatz for the scalar field 

= e-''''+'"'^Y{e)R{r) + h.c. . 

The equation for Y{6) is the standard equation for the flat space spherical harmonics plus an extra 
term that can be neglected if {a/tf <^ 1. As before, the superradiant condition ^ implies that 
{a/iy < 1/4, so we will always use this approximation (one can check its accuracy using known 
numerical results for oblate spheroidal harmonics, see e.g. [32]). Then the equation for the radial 
function R takes the form 

Adr{AdrR) + {uj\r^ + - Aavgrmuj + a^m^ - A(/iy + a^w^ + /(/ + 1)) = . (12) 

At the horizon, a non-singular solution of this equation satisfies [23] 

R = canst ■ e-''^"'"+>* as r ^ r+ , (13) 

where is the "tortoise" coordinate defined through 

dr^ = (r^ + a^)A-^dr . (14) 



2.1 Non-relativistic approximation a// <C 1 



This approximation [T81I20]. that initially was applied for superradiant scattering rather than the 
calculation of the instability rate ^6], makes use of the separation of scales between the size of 
the cloud and the black hole horizon following from relations ([9]), (11 ). The radial equation is now 



solved in two different regimes: the near and far horizon regions. In the region far from the black 
hole horizon, r ^ r^, neglecting terms suppressed by {a/l)'^, the solution takes the same form as 
the radial wave function of the Schroedinger equation with an 1/r potential, 

2 

Rfarir) = {2krye-'''U{l + 1 - 2(/ + 1), 2kr) , (15) 

Tgk 

where U is the confluent hypergeometric function of the second kind, and k is the axion momentum, 

e = fii-u'. (16) 
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In the ordinary Schroedinger equation with an 1/r potential the spectrum of frequencies u is 
determined by requiring the regularity of R at the origin. Instead, in the black hole case one has 
to impose the regularity of the field at the horizon — the incoming wave boundary condition (13). 
One way to do this is to solve the radial equation (12) in the near horizon regime. After dropping 
terms of order a/l the solution in the near-horizon region, < r — r+ ^ {l/aYrg, that satisfies 



the boundary condition (13) takes the form 



Rr, 



2Fi(-/,/ + l,l + 2?P, 



(17) 



where 



P = 2r_ 



u — raw A 



and 2F1 is the Gauss's hypergeometric function. At a/l <^ 1 the ranges of validity for the two 
approximate solutions (15) and (17) have an overlap, so the approach in [181 [20] was to match the 



lowest terms of the Taylor expansion for Rfar at small r with the asymptotic behavior of (17) at 
large r, with the following result for the imaginary part of the frequency. 



Imn 



where 



Imn 



2/ia^'+VH 



24^+2(2/ + n + l)! 
(/ + n + 1)2'+%! V(2/)!(2/ + 1)! 



[raw A 



2 I 



- /ia)C/ 



Imn 



:i8) 



+ Ar'i{mw+ - Ha)'- 



Note that the real part is well approximated by the hydrogen spectrum ([T]). These approximate 
formulaes exhibit many of the features of the full answer. In particular, the sign of Timn is 
determined by the sign of {rawj^ — f^a), in agreement with ^ within the accuracy of the non- 
relativistic approximation used to derive (18). Also, in the regime of apphcability of (18), a/l ^ 1, 
widths Timn drop exponentially as / increases. This implies that for a given value of a the fastest 
superradiant level is the one with a smallest possible I, i.e., the I = ra level with ra chosen in such 
a way that the superradiance condition ^ is satisfied. 

It will be important in what follows that as I grows the radial quantum number n for the 
fastest superradiant level grows as well. For instance, using (18) we find that for I = ra = 4 the 
instability rate for the n = 1 level is faster than for the n = level. 



440 



441 



0.9. 



(19) 



The result ( 18 ) is often referred to as the "low mass" or small a approximation. However, as 



the above derivation shows, the approximate solutions (15) and (17) hold and have an overlapping 



regime of validity at large a as well, as long as a/l ^ 1. It is true, though, that (18) may be 
quite inaccurate close to the superradiance boundary a < w+, both due to the inaccuracy of the 
approximate solutions and because the overlap interval where both of the solutions hold shrinks. 
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Figure 5: Superradiance rates obtained using our semi-analytic method (solid lines), non- 
relativistic approximation (dashed lines) and WKB approximation (dotted line) for a near- 
extremal black hole, a/vg = 0.999. Different colors correspond to superradiant levels with different 
values of the angular quantum number I. 



To have better control of the precision and to improve on the latter defficiency of ( 18 ) we adopt 
the following semi-analytic procedure, which is similar to that in |33]. Instead of matching the 
leading terms in the asymptotic expansions of (15), (17) we numerically matched these functions 
and their first derivatives at a point r^, within the overlap region. To find the optimal value for 
we calculated the relative residuals after plugging Rfar and Rnear in the original radial equation 
(12). We pick r* as the point where the two residuals are equal. We present our semianalytical 
results for the instability rates in Fig. [s] (solid lines) together with the instability rates given by 
(18) (dashed lines). Different lines correspond to different I = m levels, and we picked the radial 
quantum number n to maximize the superradiance level in each case (the dependence on n is very 
mild) . 

These results were obtained for a near extremal black hole with a/vg = 0.999. The instability 
rate decreases very slowly with spin for small a, and the main effect of reducing the spin is that 
the instability shuts down earlier, at am = mw+{a) for the different / = m levels, as follows from 
the superradiant condition (jsj) . To illustrate this in Fig. [6] we present the superradiance rates for 
/ = 1 level for several values of a/rg. 

We see that (18) perfectly agrees with our semi-analytic results at sufficiently small a/l, but 
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Figure 6: Superradiance rates obtained using our semi-analytic method (solid lines) for different 
values of the black hole spin. 



they quickly start being different indicating that the procedure of matching the leading terms of 
the asymptotic expansions that results in (18) is not very accurate. Of course, our semi-analytic 
results are also not good enough for precision calculations close to the superradiant boundary, 
a ~ mw^, however, they agree quite well with numerical calculations in [30] (at least for I = 1, 2, 3 
presented in [30] ) and with the WKB results presented in the next subsection. In particular, the 
superradiance rate of [30] is maximum, F ^ 1.5 ■ lO^^r"^, at a ~ 0.42 in a good agreement with 
our results. As we will see, many of the observational consequences of superradiance are not very 
sensitive to the exact values of the superradiance rates at a ~ mw+, so this level of precision is 
enough for our purposes. 



2.2 WKB approximation a ^ 1 

Another useful approximation for the superradiant rates [19] , complementary to the slow velocity 
expansion above, is the WKB method that can be applied at a S> 1. We will closely follow the 
methodology of [19], however, our results disagree with [19] by an important factor of two in the 
expression for the tunneling exponent. Most likely, this factor was accidentally missed in [19] (this 
discrepancy was also pointed out in [31] without any derivation). 

In this approach, the tail of the wave function that propagates in the ergo-region, where su- 
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Figure 7: The shape of the radial Schroedinger potential for the eigenvalue problem in the rotating 
black hole background. Superradiant modes are localized in a potential well region created by 
the mass "mirror" from the spatial infinity on the right, and by the centrifugal barrier from the 
ergo-region and horizon on the left. 



perradiant amplification takes place, is calculated using the WKB approximation. This is just a 



classic tunneling calculation. Indeed, after switching to the tortoise coordinate ( 14 ) and introduc- 



ing = (r^ + a'^Y^'^R the radial equation (12) takes the form of the Schrodinger equation 



with the potential 



2 Argramuj — a^m^ A f ^ l{l + I) + k'^a^ Sr^ — AvgV + a? SAr^ 



^ + (^2 + ^2)2 + r2 + a2 ^ r2 + a2 ^ {r^ + a2)2 {r^ + a^f 

(21) 

We include the (— w^) term in the definition of the potential, because even if we were to separate 
it, there would be a residual dependence on uj. We present the qualitative shape of the potential 
V for a typical choice of parameters in Fig. [7| One can clearly see the potential well where the 
bound Keplerian orbits are localized and a barrier separating this region from the near-horizon 
region where superradiant amplification takes place. 

Consequently, the axion wave function at the horizon r = r+ (corresponding to r* = — oo) 
is suppressed relative to the wave function in the vicinity of the Keplerian orbit by a tunneling 
exponent, 

|i^(r+)|^|i^(r,)|e-^ 
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where the tunnehng integral / is 



dnVv= / cir ^^^" (22) 

with ri 2 being the boundaries of the clasically forbidden region. We will only follow the leading 
exponential dependence on and do not aim at calculating the normalization prefactor in front 
of the exponent. 

To relate the tunneling exponent with the rate of superradiance instability let us consider again 
the energy flow equation ([6]). Integrating it over the horizon we obtain 



uj{mw+-io) I \Y{e)R{r+)W (23) 

'horizon 



dt 



where £ is the energy in the axion cloud. The energy is maximum in the Keplerian region, so that 
in the limit where we only keep track of the dependence on the exponent we can write 



8 oc \R{r^)\' ~ e^^\R{r 



and, consequently, to rewrite (23) as 

— = const ■ {mw+ — u)e £ . (24) 

In other words, the WKB approximation for the superradiance rate give^ 

r = 7(mw+ — w)e~^^ , (25) 

where the normalization prefactor is determined mainly by the spread of the wave function in the 
classically allowed region. We will limit ourself by calculating the exponential part F. We leave 
the technical details for the Appendix, and present only the final result here. Namely, the final 
answer for the tunneling integral in the extremal Kerr geometry takes the form 

I = n(2a - v^2a(a - 1)) , (26) 

which translates in the following superradiant rate, 

TwKB ~ 10-\-ie-2™(2-v^) ^ , (27) 



where we took the large a limit in (26) and chose the prefactor to match the low a results of 
section 23] (this value also agrees with that of [19] and [S])- As we already said, the exponent 



in (27) is larger than that in [19] by a factor of two. As explained in the Appendix, the rate (27) 
provides an upper envelope for superradiance rates at different / in the large a limit. We have 
presented (27) by a dotted line in Fig. [sj it agrees reasonably well with the previous a/l <^ I 
results. 



^Notc, that at this stage we still agree with [19]. 
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3 Dynamics of superradiance 



Let us turn now to discussing the dynamical consequences of the superradiant instabihty. One 
important property of the rates calculated in section [2] is that the time-scale for the development 
of the instability is quite slow compared to the natural dynamical scale close to the black hole 
horizon, F"^ > lO^r^. Consequently, in many cases non-linear effects, both gravitational, and due 
to axion self-interactions, become important in the regime where the system is still in the quasi- 
linear regime, so that non-linearities can be treated perturbatively. Then the dynamics of the 
axion cloud can be described by the following set of kinetic equations for the occupation numbers 
Ni for different levels, 

^ = TijNj + T^jkN.Nk + ... . (28) 
This equation gets simplified in the quasilinear regime where we can truncate the expansion in 



the r.h.s. of (28), by keeping just a finite number of terms. Note that for simplicity, we drop 



A'^j-independent terms in the r.h.s. of (28); those terms correspond to spontaneous emission. 

To avoid confusion, let us clarify the following. Throughout this paper we often use quantum 
terminology (occupation numbers, transition between levels, etc.) to describe the axionic cloud. 
This appears to be perfectly appropriate given that the size of the cloud is comparable to the 
Compton wavelength of the axion. On the other hand, occupation numbers for all dynamically 
relevant levels will always be exponentially large in what follows, ~ 10™, so that all the dynamics 
can be accurately described by a classical field theory. Of course, there is no contradiction here and 
both descriptions are correct. The very fact that we can use classical field theory to describe the 
dynamics of axion particles in the cloud reflects its quantum mechanical origin. For instance, in 
the classical field description instead of using occupation numbers Ni one can Fourier decompose 
classical field into harmonics with different frequencies and follow the (squared) amplitudes for 
different harmonics. Given numerous analogies with atomic physics we find the quantum language 
useful in many cases, but will also use the classical one, when more convenient. 



Coming back to the kinetic equations (28), at the linear level, the r.h.s. of (28) is determined 
by the superradiant rates presented in section [2j Namely, 

Tij ^ij^i ; (^9) 

where Fj is the imaginary part of the frequency for the i-th. level — positive for levels satisfyin g ([3] ) 



and negative otherwise. There could be other model dependent sources of linear terms in (28). 
For instance, if an axion has an electromagnetic coupling 

e^'^'^'F^^F^p , (30) 



where C is an order one constant (for the QCD axion C = 4/3 in 4d grand unified theories), then 
in the presence of a magnetic field axions will convert into photons with a rate [1] 
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where the reference values for the parameters are chosen to be those for the QCD axion, and the 
choice of a reference magnetic field is motivated by the estimate (Eq. (75) of [35] ) 



E ~ 4 X lO^G — (32) 
for the largest magnetic field the accretion disc can support near the horizon of a black hole of 



mass M. It is evident from (31) that axion-photon conversion is too slow to be relevant for the 
dynamics of superradiance. 



To describe the development of the superradiant instability one needs to supplement (28) with 
equations for the time evolution of the black hole mass M and spin J. These depend on the 
environment of each individual black hole, but, in general, accretion in the absence of mergers or 
other violent events has a characteristic time scale whose lower bound is set by the Eddington 
time te 

The value of te shows that the superradiant instability time, even though much shorter than the 
black hole infall time, is much faster than the evolution time scale for astrophysical black holes. As 
a result, for most of the discussion that follows we will ignore accretion, unless stated otherwise. 



For example, in section 4.5 we discuss a particular accretion model. 



Let us describe now the leading non-linear processes, which should be included in (28) to 
describe the development of the superradiant instability. The two sources of non-linearities are 
gravitational interactions of axions and non-linearities in the axion potential itself. 



3.1 Gravity wave emission 

The axion cloud may lose its energy and angular momentum by emitting gravitational waves. 
There are two major processes giving rise to graviton emission. The first is analogous to photon 
emission from atoms — this is just the transition of axions from one level to another. The major 
difference in the present case is that the transition rates between populated levels get enhanced 
by their occupation numbers. 

The other process of graviton emission is less familiar. Unlike electrons in the atom, axions 
do not carry any conserved charge. Consequently, they can emit gravitational waves also through 
processes that do not conserve the axion number. In particular, one-graviton annihilation of two 
axions is possible at the same order of perturbation theory as the transition between different levels, 
i.e., the corresponding amplitude is proportional to Mp^. Of course, this process is kinematically 
forbidden in flat space, where only two-graviton annihilation (with amplitude proportional to Mpf) 
is compatible with energy and momentum conservation. However, the black hole gravitational 
field breaks translational invariance, so that one-graviton annihilation is allowed for axions in the 
presence of a black hole. The closest analogue of this process in atomic physics is the one-photon 
annihilation of a positron with an atomic electron [36] . 

This annihilation process is quite unconventional, because the graviton momentum is deter- 
mined by the axion mass kg ^ 2iia and is parametrically larger than the momentum of axions in 
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the cloud ka ~ {a/l)fia- In other words, unhke for conventional astrophysical sources of gravita- 
tional waves, the wavelength of an emitted radiation is not parametrically longer than the size of 
the source (even though the cloud is non-relativistic), and the standard quadrupole approximation 
does not apply. 



3.1.1 General formalism 

The calculation of the gravitational wave flux due to both kind of processes is further complicated 
by the fact that the system is bound by gravity, and the axion cloud is not that far from the 
black hole horizon, where the gravitational field is non-linear. However, these complications can 
be safely ignored for order of magnitude estimates, which are enough for our purposes. To perform 
the estimates it is convenient to switch to the classical field theory description. Let us write the 
axion field of the cloud in the form 

4> = J2e-'^'x{u,^,r,e) + h.c. , (34) 

where the frequencies u ran over different bound levels in the black hole geometry. The related 
expansion of the scalar energy-momentum tensor takes the form 

T^u = J2 e~"'*^M-(^^'> ^, r, 9) + h.c. . (35) 

To estimate the gravitational wave flux from the axion cloud we use the flat space formula [37] 
for the gravitational wave power, 

dP GnOo'^ rprp J.J. 

where 

r,j(w,fc)= /"rf=^x7^j(w,x)e"*^'' (37) 



where x denotes the flat space Descartes coordinates, \k\ = u, and the TT superscript stands for 
the projector on the transverse-traceless part, 

where 

_ _ kjkj 

In this classical language the graviton emission due to axion transitions between levels correspond 
to terms in the harmonic expansion (35) of the energy- momentum tensor of the form 
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where, for simplicity, we dropped all the derivatives appearing in the expression for T^y. On the 
other hand, the annihilation processes correspond to terms of the form 

e-<-+-'^'x{uj)x{uj') + h.c. 

To calculate the gravitational wave power, one needs to solve for the scalar field harmonics x{^) 



in the black hole background and plug them into (36). This is similar to calculating the trajectory 



of a compact stellar mass object falling into a galactic mass black hole in the probe particle 
approximation, and using this trajectory as a source in the linearized Einstein equations to solve 
for the hij components of the metric, that determine the flux of gravitational waves. In principle, 
the latter step should be done in the curved geometry — the flat space approximation is just an 
estimate both for the overall rate and for the frequency profile of the gravitational wave signal. 
For instance, for the annihilation signal from a single populated level with the frequency u the fiat 



space expression (36) predicts the monochromatic gravitational wave line of frequency 2a;, while 
taking into account deviations from the fiat space would induce the broadening of the line due to 
gravitational redshift. 

However, given the separation of scales between the size of the axion cloud and the black hole 



size, we expect that by using the fiat space expression (36 ) we are making at most order one mistake 
in the overall rate, which is accurate enough for our purposes. For the same reason, the spectral 
distortions should be rather small. It is likely that one still needs to calculate more accurately 
the spectral shape of the signal (the waveform) to be able to observe it with gravitational wave 
detectors. This is a technically involved calculation, which is beyond the scope of the current 
paper. Note that in [2H], where gravitational waveforms in the Kerr metric are calculated, the flat 
space approximation gives quite accurate results. 



The other approximation we adopt to evaluate (36 ) is rather than solving for the eigenfunction 
x(w,<^,r, ^) in the full Kerr background we will just use the Newtonian approximation for the 
metric. Again, this is justified because the cloud is mainly localized relatively far from the black 
hole horizon. In this approximation the eigenfunctions xi'^i ^i ^i ^) ^^^e the familiar wave functions 



for the electron states in a hydrogen atom, and the evaluation of (36) becomes straightforward. 



Let us present here some representative results of this calculation that we will use later in section |4} 
3.1.2 Axion transitions between levels in a black hole atom 

Let us start with the more familiar case of graviton emission due to axion transitions between 
levels. As we will see in section |4] an important source of potentially observable gravitational 
radiation is related to axion transitions between levels with equal angular quantum numbers / and 
m, but with different radial quantum numbers n. This process will be relevant for sufficiently 



high Z, such that the fastest superradiant level has n > 0. For instance, as we saw in section 2.1 



for I = m = A level the n = 1 superradiant level is faster than n = 0, Eq.(19). For the transition 



between these two levels (36) gives 



^^(6, 5g) ^ iViiVo '''^'\t^r" sin^ ^ + ■ ■ ■ ^ 3 • lO-^^AT.ATo^ sin^ 9 , (3 
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where A^^i and Nq are the occupation numbers for the two levels and dots stand for higher order 
terms in the small a expansion (we checked that they can be neglected in the superradiant regime 
«<2). 

For transitions between levels the wavelength of an emitted graviton is much longer than the 
size of the system, so the conventional quadrupole formulae for a graviton emission should also 
be a good approximation. As a cross-check of our calculation let us compare result (38) with the 
quadrupole formula. For the transition rate between two levels the latter gives [3^ 



dNi 



dt 



quadr 



where Ao; is the frequency splitting, which in our case is equal to 

-2/1 1 



25 36 



and lij is the transition mass quadrupole moment. For the estimate we take 



2 4 



where the size of the cloud Tc is estimated by (11). Altogether, this gives 



dNi 



dt 



{<6g ^ 5^) ~ 8 ■ 10" WiA/'o- 



quadr 



On the other hand, from (38) we get 



dt Aw 



^ 3 ■ lO-'iViiVo^V 



(39) 



(40) 



(41) 



(42) 



in a perfect agreement with (41). 

Proceeding as above it is straightforward to calculate other transition rates. For instance, 
for transitions from the fastest (ra = 0) / = 2 superradiant level to the fastest (also n = 0) 
I = 1 superradiant level one gets (for simplicity, we present only the total power integrated over 
directions) 

P(3i ^ 2p) « iV,7V„^I^^^^ « 4 . 10-7V,/V„^ . (43) 

This rate is suppressed by a higher power of a, because the quadrupole transition between these 
two levels is forbidden. 



3.1.3 Axion annihilations 

We see that for estimating the transition rates one can use the standard multipole formulae. As 
we said, this is no longer the case for annihilations, where the wavelength of an emitted graviton 
is shorter than the size of the cloud. In this case the suppression for the emission rate is related to 
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Figure 8: Total annihilation rates in units of r„ ^ for ^f^^^ = 10 ^ as a function of all for the 

^ 9 Mbh ' 

superradiant levels with / = 1, 2, 3. 



the decoupling of high momentum modes — the Fourier transform in (37) involves convolution of 
the slowly varying energy-momentum tensor of the axion cloud with a rapidly oscillating exponent. 
A direct calculation gives the following result for the annihilation rate at the 2p level, the fastest 
superradiant level, 
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(2 X 2p ^ graviton) ^ A^^ ^^^i^" (35 + 28 cos 29 + cos AO) , (44) 



where we again expanded the full answer at small a and is the occupation number. Note 
that, unlike for transitions, the corrections from higher order terms in a change the answer for 
the annihilation rate by a factor of order one close to the upper boundary of the superradiant 
regime ([3]) (a ~ 1/2), however this level of precision is enough for our estimates. We see that the 
suppression for the annihilation rate at small a is much stronger than for transitions. 

For higher I levels the suppression at small a becomes even stronger, because the multipole 
number of the emitted graviton grows with /. For instance, the annihilation rate for the 3(i level 
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scales as a"^^. Still, given that the power of a is already very high even for the p-level annihilation 
and the suppression in the annihilation rate is determined by how non-relativistic axions in the 
cloud are, one may expect the annihilation rates for different levels to be comparable at equal 
values of the axion velocity, or, equivalently, at equal values of a/l. To illustrate that this is 
indeed the case, we present in Fig. [8] the annihilation rates for the first three superradiant levels, 
/ = 1, 2, 3, in units of Vg for a cloud mass of 10"*^ Mbiack hoic as a function of a/l. We see that they 
indeed agree within an order of magnitude. 



3.2 Axion non-linearities 



Let us turn now to another important source of nonlinear terms in (28) — self-interactions of the 
axion field itself. As we discussed, the superradiant cloud is always close to be non-relativistic, 
so let us first discuss the self-interaction effects in the non-relativistic approximation. In the 
non-relativistic limit, the axion field takes the form 

(e-^'^-V + e^'^^V*) , (45) 



where the characteristic scales for space and time variations of the function are much longer 



than /x^ ^. Then we plug the ansatz (45) in the axion action and drop all rapidly oscillating terms. 



As the result we obtain the following effective action for tjj, 

Snr = j d^x [irdt'il^ - ^d.ijd.r - + , (46) 

where $ is the Newtonian gravitational potential of the black hole, and we kept only the leading 



non-linear term from the axion potential. If we also drop the quartic self-interaction term in (46) 



then the equation following from (46) is the conventional Schroedinger equation in the external 
gravitational field. As before, the superradiant instability can be thought of as coming from 
an unconventional boundary condition at the origin. With the quartic term taken into account. 



the axion action (46) leads to the non-linear Gross-Pitaevskii equation, well-known in condensed 



matter physics to describe the dynamics of an interacting BEG (see, e.g., [39j for an introduction). 



The sign of the interaction term in (46) corresponds to an attractive interaction between axions. 
The black hole gravitational potential plays the role of the BEG trap. 

Note, that as compared to gravity, axion self-interactions give rise to higher-order non-linearities 
(quartic, rather than cubic), however, they are suppressed by the scale fa, which can be signifi- 
cantly lower than Mpi. 

Let us describe the major consequences of these non-linearities. We postpone the systematic 
discussion of how superradiant instability develops till the next section |4| but it is intuitively clear 
that typically the axion field in the cloud is dominated by a single harmonic in the expansion 



(34) — the one corresponding to the fastest available superradiant level. This is especially natural 
to expect at small and moderately large values of a, when superradiant levels are very sparse. So 
let's first consider how non-linearities affect a cloud composed of a single superradiant level. 
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3.2.1 Bosenova 



As the number of axions in the cloud increases, the attractive force between axions becomes more 
and more important and at some point the shape of the cloud changes significantly as compared 
to the one corresponding to the unperturbed hydrogen wave functions. We can estimate when 
this happens by equating the potential energy of axions in the cloud to the self-interaction energy, 

^ ~ ^ . (47) 



By integrating (47) over the volume we obtain that self-interaction effects become important in 



determining the shape of the cloud when 



where N is the number of axions. Here we made use of (11) for the size of the cloud, and set 



n ~ /. It is more convenient to write this bound as a condition on the mass of the axion cloud. 

We see that for typical values of the parameters we are interested in, fa ~ Mqut, non-linearities 
start playing an important role in determining the shape of the cloud quite early — when the cloud 
constitutes only 10~^ 10~^ of the black hole mass, and even earlier if the axion scale fa is 
significantly below the GUT scale. 



As the cloud grows and its size is close to saturating (48), the shape of the cloud is deformed 
and is no longer determined by the hydrogen wave functions. However, a much more dramatic 
effect happens as the size of the cloud keeps growing. The effect was experimentally observed 
in trapped EEC's with attractive interactions and is known under the name "Bosenova". Above 
some critical mass of the cloud the gradient energy of the axion field ( "quantum pressure" ) cannot 
compete with the attractive force due to self-interactions and the cloud collapses. Indeed, it follows 



from (46) that the energy of the static cloud has the following parametric dependence on its size 



V(R) ^ N'-^^ - N- + ^ . (49) 



At small this energy has a minimum corresponding to a (meta) stable cloud, however, at large 



A^ the last term in (49) dominates over the repulsion due to the quantum pressure term, and the 
cloud collapses. 



3.2.2 Shutdown of superradiance due to level mixing 

Yet another important consequence of non-linearities is that they may stop further development 
of superradiance. To see how, let us first consider the Press-Teukolsky "black hole bomb" — a 
rotating black hole surrounded by a spherical mirror. As we discussed in the Introduction, this 
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system provides the simplest example of a superradiant instability — a single photon introduced 
inside the mirror and satisfying the superradiance condition ^ bounces between the mirror and 
the horizon and gets amplified in the ergo-region. Now imagine that the mirror has a defect and 
does not possess a perfect rotational symmetry around the axis of black hole rotation. Then, when 
scattering off the mirror, some of the photons satisfying the superradiant condition ([s]) change 
their quantum numbers. As a result they may exit the superradiant regime and be absorbed at 
the horizon. If the defect is substantial enough this may dump the superradiant instability. 

Coming back to axion superradiance, as a consequence of self- interactions the axion cloud itself 
acts as a defect and may dump the further development of the instability when it becomes large 
enough. To analyze this at a more quantitative level, let us write the axion wave function in the 
form 

tp = ^o + 5'(p , 

where ^jJo oc g-^C'^ot-mov) jg i^j^e field of the axions populating the most occupied level, and Sip 
is a perturbation. We assume that the time scale for the growth of the cloud is much longer 
than the oscillation period for ijj and consider the dynamics at scales short compared to the 
instability time — under these assumptions uo can be taken real. Then the linearized equation for 
the perturbation 6ijj has the following form, 

d'^Sib 1 
idtSiP = -^+ IJ^a^ - ^ {2roi^oSip + ijl^r) . (50) 

The interesting feature of this equation is that it mixes ip and ip*; this property is typical for BEC 
perturbations and gives rise to the notion of Bogoliubov's quasiparticles. Namely, the solution to 



(50) mixes postive and negative frequency components, 

Sip = e-*('^o*-"^o^) (u(r, ^)e-*('5'^*-'5™^) - v*{r, e)e''-^^^-^'^^^) . (51) 

In the next section we will be interested in a situation, when the unperturbed cloud has parameters 
very close to the boundary of the superradiant region, uq = uiqw^ and the perturbation 5%p 
corresponds to the fastest available superradiant level, which is / = m = mo + 1. This same 



reasoning applies for the levels with m higher than mo + 1. The first term in (51) would give rise 
to such a perturbation, however, we see that as a result of the interaction with the background 
BEC, the perturbation has also an admixture of the non-sup err adiant m = mo — 1 modes. Of such 
modes the one with the fastest damping rate also has / = mo — 1, and as far as we can tell, there 
is no reason that would forbid an order one overlap of the function v with the I = m = mo — 1 
mode. 



To see whether the perturbation (51) is superradiant or dumped, let us proceed as we did 
before in the derivation of the time averaged energy flux through the black hole horizon, (|6]). The 
axion field has now three different harmonics with frequencies Uq, Uq ± 6u. However, only the 
latter two contribute to the flux, because the first one saturates the superradiance condition (and 
all cross-terms vanish as a result of time-averaging). As a result, the flux takes the following form 

{P^lG^) = wi(u;i - (mo + l)w+)\uh\^ + ^2(^2 - (mo - l)w+)\vh\^ , (52) 
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where uji^2 = /^a + ± Su, and Uh and Vh are the values of the functions u and v at the black 



hole horizon. For our choice of parameters, the u-term in (52) gives rise to the energy flux from 



the black hole, while the f-term gives rise to the flux into the black hole. To deduce the direction 



of the net energy flux, let us recall that derivation of the WKB formula (25) implies, that the 



functions Uh,Vh at the horizon are related to their values Uc,Vc at the location of the cloud as. 



\Uh\ 
\Vh\ 



£1 

r2 



Ur 



where Fi and F2 are superradiance and dumping rates for the two levels. In turn, the ratio 
Uc/vc is determined by the relative strength of the non-holomorphic in Sip term in the perturbed 



Gross-Pitayevskii equation (50) 
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Combining these two relations together, we obtain that the perturbation Sip is superradiant if the 
number of axions in the cloud is smaller than 



£1 

r2 



1/2 
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As in ( 48 ) it is convenient to rewrite this condition as a bound on the fractional mass of the axion 
cloud, 



BH 



< 



1/2 



)_J_a_ 



(53) 



To estimate the ratio of rates in (53) we can use (18). In Fig. 9 we plot the result for the first 



few values of mo as a function of a/vg. We see, that for a broad range of a/vg this ratio changes 
between ~ 10~^^ and ~ 10"^. Note that for a/rg close to one the non-relativistic approximation is 
not accurate for the / = itlq — 1 level. Comparing with the numerical results of |30j for mo = 1,2 
suggests that the ratio IF1/F2I is actually close to 10~^ at a/vg close to one, rather than to 10^^ 
as shown in Fig. [9j 



By comparing (53) and (48), we conclude that a large occupation number for one of the 



superradiant levels may indeed strongly inhibit the development of superradiance for other levels 
even in the regime when non-linearities are still irrelevant for the shape of the cloud. 

Importantly, this does not happen for the most occupied level itself — this corresponds to 
considering 5w = 5m = in (51). The derivation of the horizon flux (|6|) did not assume that 



the field is linear, so that as soon as the field is well approximated by a single exponent ipo oc 
^-i(u)ot-mo<fi) ^Yie superradiance will continue even in the non-linear regime. 

It is important for this argument that the ansatz i/jq consistent solution 

of the Gross-Pit aevskii equation at the non-linear level. This is no longer true if one considers 
the full scalar equation including all relativistic corrections — the higher harmonics get generated. 



However, we do not expect those to change the conclusion. As shown in section 3T, relativistic 
processes are strongly suppressed for superradiant levels. Still, this point deserves further study. 
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Figure 9: Ratios of the superradiance rate for the (Z + 1) level to the absorption rate for (/ — 1) 
level if the / level has a vanishing imaginary part for / = 1, 2, 3. 



Note, that the effects discussed so far — deformation of the shape of the cloud and shutdown 
of superradiance due to level mixing — in principle could be caused by gravitational backreaction 



of the cloud as well. We did not discuss them in section |3.1[ because the shape deformation is 
always small as soon as Ma <C Mbh, and the level mixing is absent in the non-relativistic limit 
when the density of axions pa ~ fJ'a4'*4' does not depend on the azimuthal angle ip, as the cloud 
is dominated by a single level. Consequently, these gravitational effects are likely to be always 
subdominant with respect to those caused by axion self-interactions. 



Conversely, axion self-interactions may cause effects similar to those discussed in section 3.1 
annihilations and elastic scatterings of axions. The leading annihilation process is annihilation 
of three axions from the cloud into one axion in the continuum. Given that the outgoing axion 
has energy of order 3/ia, in the leading approximation it can be considered as massless, and the 



calculation of the emission rate can be done similar to the graviton case. The analogue of (36) for 
three axion annihilation in the massless approximation reads 



dP 



sin OdOd^p 

where 



2tu^\jicu,k)\\ (54) 
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Then a straightforward calculation gives the following result for the annihilation of three axions 
from the I = m = 1 level into a single unbound axion, 



dP 2^°a23 

-(3 X Ip — )■ continuum) ~ iV^ ^— ^ sin^ 9 . (55) 



To compare the efficiency of this process to the one graviton annihilation (44), let us integrate 



over the angles in both cases to calculate the total emitted power and take the ratio. We get 

P(3 X Ip — >■ continuum) ^^_2 4 Mp^ 
P{2 xlp^ graviton) " Mbh^ ' 

so that self-interactions dominate when the size of the cloud is not too small, however, as the cloud 
decreases, the two axion annihilation into gravitons takes over. For instance, it will typically be 



more important when the cloud size approaches the value in (53). 

Just as in the case of graviton emission one may look for the processes that conserve axion 
number and as a result may be less suppressed by powers of the small axion velocity in the cloud. 
An obvious candidate process is an elastic two-to-two scattering of axions. However, it appears 
likely that such processes are not important. The reason is that for the two-to-two scattering to 
be unsuppressed, three of the participating axions should correspond to highly occupied levels. 
As we will discuss in section |4| typically at any given moment of time only few of the levels are 
significantly occupied, and usually there is one which dominates the cloud. Then the most likely 
candidate for the scattering process is a scattering when two axions from the most populated level 
scatter, and one goes down to another highly populated level, while the remaining axion files out 
in the continuum. For this scattering to be compatible with the energy conservation one needs 

2 1 

ni n^ 

where ni and n2 are the principal quantum numbers of the high- and low-lying levels, respectively. 
The discussion of the superradiance development in the next section implies that the situation 
where two levels satisfying this condition simultaneously have large occupation numbers is hardly 
possible. 

Note that throughout most of the discussion of axion self-interactions and gravitational wave 
emission we treated axions in the cloud as free particles with hydrogen wave-functions, while 



strictly speaking the elementary excitations of the axion BEC are Bogoliubov quasiparticles (51). 
The free particle approximation is accurate when the mass of the cloud is small, so that the 
effects of self-interactions are weak, but may be misleading when the mass of the cloud is close 



to saturating the bound (48). Unfortunately, it's hard to improve on this approximation without 
going into numerical simulations of the cloud, which are beyond the scope of our paper. We 
proceed under assumption that the free particle approximation is a reliable guide for an order 



of magnitude estimates at the masses close to saturating the bound (48) as well. The level 
mixing phenomenon described above provides an important example of a situation, when the free 
particle approximation is not adequate even at very small masses of the cloud. There is a clear 
physical reason why this happens — some of the levels have damping rates orders of magnitude 
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faster than the population rates for the relevant superradiant levels, so even a tiny level mixing 
qualitatively changes the dynamics. Fortunately, this effect is straightforward to take into account 
perturbatively, as we did. 

To summarize, this discussion implies that self-interaction effects cause a strong influence on 
the phenomenology of superradiance. The proper taking into account of these effects is one of the 
major challenges for obtaining an accurate quantitative description of superradiance development. 
The estimates presented here are far from being a complete accurate treatment and it appears 
likely that numerical simulations are required to really solve the problem. It is worth mentioning 
that level mixing can also be caused by the accretion disk of the black hole or a massive object 
orbiting the black hole but we have already shown in [Ij that these can be safely ignored. We 
proceed now with estimating the possible observational signals of superradiance. 



4 Observational signatures 

Now that we are well equipped with the details of how superradiance works, let us put them 
together and develop an intuition about the way superradiance develops in realistic environments 
and about the observational signatures we may expect from this process under various circum- 



stances. The full treatment of the set of kinetic equations describing superradiance (28), (70), (72) 
appears to be quite challenging, given that a large number of competing processes with drastically 
varying time scales is involved. Our strategy will be to start with a highly idealized situation 
including a minimal number of dynamical ingredients and then keep adding more processes to get 
closer to a realistic description. We already presented the list of possible observational signatures 
of superradiance in the Introduction. Clearly, following the above strategy the very first signa- 
ture to discuss is the absence of rapidly rotating black holes of size matching the axion Compton 
wavelength — the black hole spin-down is the most direct consequence of superradiance. 

4.1 Black hole Regge trajectories 

To get a rough idea of the bound on the axion mass that could come from black hole spin 
measurements, we present in Fig. |2] regions in the black hole mass vs axion mass plane where the 
superradiance timescale for a maximally rotating black hole is shorter than the age of the Universe 
and the Eddington accretion time. We used the superradiance rates of section |2] to produce this 
plot. For superradiance to have a non-negligible effect on the black hole spin the process should 
last for many e-foldings for the produced axions to carry away a noticeable fraction of the black 
hole spin. To estimate the required number of e-foldings, note that, approximately. 



Ha 01 \ Mq 



of axions need to be produced for their total mass (spin) to be of order the black hole mass (spin). 
This requires ~ 10^ e-foldings of superradiance; we took this factor into account in Fig. |2| by 
presenting the region where the age of the Universe (or Eddington accretion) time is longer than 
hundred superradiance times. We see from the plot that, as the black hole mass grows, the size 
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of the interval in axion mass which the black hole can in principle probe shrinks, because the 
superradiance time gets longer. 

This plot is useful as a zeroth order estimate, but cannot be used to deduce limits on axions 
from data on existing black holes, since the black holes being observed do not all have spins close 
to the maximum. Instead, we need to know the regions in the black hole "Regge" plot, the spin vs 
mass plane for black holes, where we do not expect to find black holes if an axion with a certain 
mass exists, as shown in Fig. |3j The intricate structure of these plots is a manifestation of level 
quantization in the superradiant gravitational atom, and their underlying physics is explained in 
what follows. 

To start with, let us ignore accretion and consider a black hole that starts off as maximally 
rotating. This approximation should be physically relevant for stellar mass black holes produced 
as a result of a fast catastrophic event, such as the supernovae explosion. Deviations from spherical 
and axial symmetry are believed to be crucial for supernovae explosions, so there should be a lot 
of angular momentum available when the black hole forms, and one may expect high initial values 
for the black hole spin — as soon as the supernovae core gets rid of all the angular momentum 
above the extremal value it collapses and forms a rapidly rotating black hole. This expectation 
seems to be supported by observations — for instance, the high value of the spin-to-mass ratio 
^ ^ 0.92;!;qq7 deduced |1D] for the black hole primary in the extragalactic X-ray binary LMC X-1 

is hard to reconcile with the young age (~ 5 ■ 10^ yr) of the system, if the spin were not natal. 

For concreteness, we assume that the parameter a for the black hole is small, a < 1/2, so 
initially the fastest superradiant level is the 2p level with I = m = 1. Then, initially one can 
ignore all the levels apart from the 2p superradiant level. It is straightforward to generalize all the 
discussion below to smaller initial values of the black hole spin and larger values of a. The black 
hole will start to lose its spin rapidly by populating this level. The time scale for this process can 
be really fast — from section 2^ we know that the superradiance rate can be as fast as lO^r^^, which 
corresponds to 10^ seconds for a two solar mass black hole. However, there is a critical value of the 
black hole spin ai(a)/rg, at which the superradiant condition ^ ceases to hold for the I = m = 1 
level, so the width of this level becomes zero and the spin-down process terminates. 

Let us for a moment consider the case where the axion self-interactions are absent, as if we 
were dealing with a free massive field rather than a (pseudo)Nambu-Goldstone boson. Then at 
this point superradiance would continue by populating the second superradiant level I = m = 2 
{3d) at a much longer time scale. Note, however, that for many e-foldings of superradiance the 
black hole spin would remain approximately constant and equal to ai{a)/rg. Indeed, if the spin 
would significantly drop below this value the frequency of the 2p-level would acquire a negative 
imaginary part, so that the black hole would start absorbing axions from the 2p-level and spinning 
up back with a rate much faster than the population rate for the 3d level. Instead, the spin stays 
practically constant close to ai{a)/rg as the black hole populates the 3(i- level while being slowly 
fed by axions from the 2p-level. When the occupation number for the latter level A^2p drops below 
^r^N^d the spin-up rate due to 2»-level cannot compete with the spin-down rate due to 2d-\eve\ 
and the black hole spin further drops down till the value 02 (a) /r^, where the superradiance rate 
for the Bd-level turns zero and the story repeats this time involving the I = m = 3 (4/) level. 

Of course, from section [3] we know that this story cannot be an accurate description of what 
actually happens — non-linearities related to the axion self-interactions and due to gravity cannot 
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be neglected in a realistic description of superradiance. However, the above simplified example 
correctly captures the major important feature — during the spin-down black hole spin tends to 
evolve rapidly till it reaches the line ai{a)/rg in the Regge plane where one of the superradiant 
levels changes the sign ("Regge trajectory"), where it can get stuck for a quite long period of time. 
In fact, as we will see momentarily, non-linear effects make this behavior even more pronounced. 

Indeed, as we discussed in section 3.2[ even relatively small amount of axions populating one 
superradiant level may shut down the instability for the next level. For instance, in the above 
example, when the black hole reaches the first Regge trajectory ai{a)/rg the 3d level does not 
start being populated until a large enough number of axions dissipate from the 2p level, so that 



its mass drops below the bound in (53). 



We discussed two processes that reduce the number of axions in the superradiant cloud, anni- 
hilations into gravitons and annihilations into unbound axions due to self-interactions. The latter 
process is more efficient at large occupation numbers. However, the annihilation rate drops down 
as the number of axions decreases and, when the cloud mass approaches the bound (53), the 



graviton annihilation, which involves only two axions, wins, as seen from (56). The total duration 



of the annihilation period before superradiance restarts, is dominated by the latest stages of the 
process. We can estimate the duration of this period by using the annihilation rates calculated 



in section 3.1 The occupation number of axions dissipating from the cloud through annihilations 
into gravitons satisfies the following equation 



dN 
~dt 



(2 X axion — )■ graviton) = —TannN 



(58) 



where the coefficient Tann can be deduced by integrating the annihilation rates, such as (44), over 



the angles. By solving (58) we obtain that the occupation number evolves in time as 
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(59) 



where at the last step we took the late time asymptotics. By using (53) we find that the anni- 



hilation time needed to clean the system before superradiance can continue to populate the next 
level is 
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(60) 



where the last step is merely the definition of the normalized annihilation time T{a). In Fig. 10 



we 



present T{a) for the first three levels. We see that for stellar mass black holes, depending on the 
parameters, there is enough time for annihilations to complete on one or two Regge trajectories. 

The above discussion gives rise to the following picture of the black hole evolution under the 
infiuence of superradiance. An initially fastly rotating black rapidly loses its spin and approaches 
the closest Regge trajectory ai{a)/rg. Then, for a long time the black hole stays at this trajectory 
while axions in the cloud keep annihilating. When the cloud mass drops below the critical value 



(53) the superradiance becomes operative again and the black hole rapidly travels to the next 



Regge trajectory. Consequently, in the region of the Regge plane affected by superradiance masses 
and spins of observed black holes should follow quantized trajectories. 
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Figure 10: The time required for the axion cloud around a 2Mq black hole to dissipate such that 
the next superradiant level can start being populated for clouds with / = 1, 2, 3,4. 



During transitions between Regge trajectories another non-linear process discussed in sec- 
tion |3.2| — bosenova — becomes important. As the cloud mass during the transition grows above 
(48) the cloud becomes unstable as a result of the attractive axion self-interactions and collapses. 
The detailed analysis of this process requires numerical simulations, which are beyond our goals 
in this paper. However, the most likely outcome seems to be that order one fraction of the cloud 
gets absorbed by the black hole and order one becomes relativistic and escapes at the time scale 
set by the size of the black hole. Condition ^ implies that to complete the transition to the next 
Regge trajectory a black hole needs to release up to a few percent of its mass into axions. Conse- 
quently, each transition proceeds through a sequence of tens to hundreds of Bosenova explosions 
for fa ~ Mgut- As discussed later, these explosions may give rise to the observable gravitational 
wave signal for supermassive black holes and perhaps to the gamma ray signal for the QCD axion. 

It is straightforward now to add accretion into this picture, at least at the qualitative level. 
Under the influence of accretion, the black hole mass and spin will still stay on the Regge trajectory, 
since the superradiance rate for the corresponding level is much faster than the accretion rate 
away from the trajectory. Indeed, if accretion brings the black hole above the Regge trajectory, 
the level acquires positive imaginary part and spins the black hole down back onto the trajectory. 
Conversely, if as a result of accretion the black hole deviates below the Regge trajectory, the 
imaginary part becomes negative and the black hole starts absorbing axions from the cloud to 
return on the trajectory. 
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Note that accretion may affect the black hole transition rate between different trajectories if 
its rate is faster or comparable to the annihilation rate. If it consistently pushes black hole above 
the trajectory, new axions will be coming to the cloud compensating the effect of annihilations. 
Conversely, by pushing the black hole below the trajectory, accretion may accelerate the dissipation 
of the cloud. 

Of course the above discussion only applies if accretion is slow compared to superradiance in 
the vicinity of the Regge trajectory. We illustrated all of the above in Fig. [3j Here, lines of different 
colors correspond to different levels. Parts of these lines where the spin increases with the mass are 
boundaries of the superradiant region, where a = mw^{a/rg) for the corresponding values of m. 
These are the Regge trajectories ai{a)/rg discussed above. For each trajectory the superradiance 
time grows at small a, and at some point becomes longer than the age of the Universe. Then, 
instead of showing the line where the width of the level is zero, we show the curve where the 
superradiance time is equal to the age of the Universe. These are the parts of the lines in Fig. |3j 
where the spin is a decreasing function of the mass. Below the solid line the superradiance time 
is longer than the age of the Universe for all unstable levels. 

To finish the discussion of the Regge trajectories, note that we started with considering an 
example of the black hole spin-down, which can be relevant for stellar mass black holes, but by 
now it is clear that also the evolution of the galactic black hole will follow the same rule — as the 
black hole enters in the region of the Regge plane affected by superradiance it starts moving there 
along the Regge trajectories, occasionally experiencing rapid transitions between different Regge 
trajectories. 



4.2 Gravitational waves 

An even more direct possibility to detect the presence of an axion cloud around black holes is 



to observe the associated gravitational wave signal. As we discussed in section [3^1] there are two 
principal processes giving rise to gravitational waves from the cloud — axion transitions between 
levels and axion annihilations. The bosenova collapse may also give rise to a burst of gravitational 
waves. 

Let us start with the transition signal. For transitions to be efficient one needs large occupation 
numbers for two different levels to be present simultaneously to get a Bose enhancement of the 



signal. Combined with non-linear effects discussed in section |3.2 this practically singles out the 



type of transitions having chances to be observed and the corresponding stages of the black hole 
evolutioi]0 



Indeed, in section 3.2 we found that even a relatively small amount of axions on the most 
populated level shuts down the superradiance for levels with different magnetic angular numbers 
as a result of level mixing. Consequently, the only chance for two levels to grow together, and as a 
result to acquire large occupation numbers simultaneously, is when the levels have equal angular 



numbers m. This case corresponds to setting Sm = in (51); the mixing for such levels does not 



change their magnetic number and does not shut down the superradiance. 



^In particular, transitions between superradiant and non-superradiant levels, chosen as an illustrative example in 
[1] , are actually never important because the transition rate is always suppressed compared to the axion absorbtion 
rate for non-superradiant levels and they never have a chance to acquire a large occupation number. 
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Furthermore, for two levels with equal magnetic numbers m but different orbital momenta I 
the superradiance rate for the more energetic level (the one with a higher I) is very much slower, so 
that by the moment the black hole reaches the corresponding Regge trajectory by populating the 
lower level, the occupation number for the higher level is tiny, and there is no significant transition 
signal. 

All this lead us to consider the transition between two levels with different principal quantum 
numbers n, but equal / and m as the most promising source of an observable gravitational wave 
signal. It is natural to consider the case, when the level with the larger principal number n has the 
faster superradiance rate. As mentioned in section 2.1[ the lowest / when such a situation takes 



place is Z = 4, so let's pick this level as the simplest representative example. The corresponding 



transition rate is given by (42). The amplitude of the gravitational wave signal at the detector is 



related to the total power emitted as 



where u is the frequency of emitted gravitons and r is the distance to the source. By making use 



of the rate (38) and plugging in the transition frequency (40) we obtain 



n ^ 10-V(e..„)V^ {^) (^) . (62) 

where ei^o are total masses of axions populating the n = 1 and n = levels, in units of the black 
hole mass, 

and the frequency of the signal v is determined by the axion and black hole masses as 

V ^ 100 Hz — ^ . 63 
\MbhJ ^ ' 

For annihilations one needs fewer conditions to be satisfied to get a significant signal — the 



occupation number for only one of the levels has to be large. As we discussed in section 3.1 the 



annihilation rates for different levels are rather similar at same values of a//, so let us consider 



the I = 1 level as a representative example. The annihilation rate (44) gives rise to a gravitational 
wave signal of strength equal to 

22 7 /^10Mpc\ fMBH\ 

where, as before, e is the fraction of the black hole mass accumulated in the cloud. The frequency 
u for this signal is given by 

u ^ 30kHz a | ) . (65) 



The numbers in (62) and (64) definitely appear interesting both at high frequencies probed by 



Advanced LIGO and corresponding to stellar mass black holes and when scaled down to low LISA 
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frequencies, corresponding to supermassive galactic black holes. However, to judge the chances 



to observe these signals we need to estimate the fractional mass of the cloud, e, in (62) and (64), 
that will determine the possible observational reach in terms of the distance to the source. 

To make these estimates note that, as we discussed, a large fraction of its time the black hole 
spends on the Regge trajectories with a relatively small axion cloud around, down to e ~ 10^^^^^", 
waiting for the cloud to be dissipated so that the next level can start being populated triggering 
a relatively fast transition to the next Regge trajectory. Transitions between different trajectories 
and relatively short time intervals afterwards, when e can be significantly larger, provide the most 
promising periods for an observable gravitational wave signal. Every black hole may experience 
several such transitions — one directly after the black hole formation, another one or two after 
periods of axion annihilations and possibly more triggered by accretion or merger events. 

From ^ we can estimate the total spin, and, as a consequence, the mass extracted from the 
black hole during such transitions. The latter typically turns out to be around a percent of the 
black hole mass. However, the main limiting factor for the size of the cloud is the Bosenova 
instability that prevents the cloud mass to grow above 10"'' -r- 10^^ of the black hole mass. This is 



enough to estimate the strength of the annihilation signal. In Fig. 11 we present the contour plot 
for the strength of the signal in the black axion mass vs a plane from axion annihilations in the 
/ = 1 level, assuming the size of the cloud e = 10"'^. We pick 20 Mpc as the distance to the source 
(which is the distance to the Virgo supercluster, containing about 2000 galaxies), and choose 
10^ seconds as a coherent integration time for the signal. Note that existing spin measurements 
(Fig. [3]) suggest a lower bound of /Xq > 3 ■ 10~^^ eV for axions matching the size of stellar mass 
black holes (which would correspond to a > 0.9 for 2Mq black hole), pushing the annihilation 
signal for Advanced LIGO into a range of somewhat uncomfortably large frequencies. 

Predicting the characteristic features of the transition signal is more involved and requires a 
detailed quantitative analysis of the dynamics, but some simple estimates can still be done. Let 
us focus on the simplest case of the 6g — )■ 5g transition. To describe the evolution of the cloud 



during the transition to the / = 4 Regge trajectory, let us truncate the system (28) by keeping 
only 6g and 5g levels and ignore accretion. Then we obtain the following pair of equations for the 
relative sizes of the 6g and 5g components of the cloud, 

^ = r44oeo - Fteieo (66) 

—7- = r44iei + r^eieo , 
at 

where r44o, r44i are the superradiance rates, and the transition coefficient Tt is determined from 
(42) to be equal to 

r, ^ 3 ■ 10-^ 



We neglected the annihilation processes which are slow compared to superradiance and transitions. 
As follows from Fig. [sjthe superradiance rates r44o, r44i are of order 10~^°r~^ for a ~ 1. Using 



the small a approximation (18) as a guide, their ratio can be estimated as r44o/r44i ~ 0.9. 



Let us focus on the case of a ~ 1. Then the dynamics following from equations (66) is quite 



simple. Both levels start being populated but the lower one has a smaller superradiance rate, and 
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Figure 11: The contour plot of constant gravitational wave signal from axion annihilations in 
the 2p level for a black hole located at 20 Mpc away from the Earth. The projected sensitivity 
curves assume 10^ seconds of a coherent integration time for LISA [H], AGIS, Advanced LIGO 
and Einstein Telescope. 



as a consequence is less occupied. By the time the occupation number of the 6g level reaches its 
maximum, ei ~ 10^^^^, the occupation of the lower 5g level is given by 



where A'^e the number of e-foldings of superradiance required to populate the 6g level. The number 



of e-foldings depends on the initial number of axions. As follows from (57) it varies from Ne ~ 165 
if initially the 5g level is not occupied, down to Ne ~ 100 if we estimate the initial occupation 
number to be determined by the dark matter density. In fact, the initial axion number can be 
significantly larger, if we consider a transition after a recent Bosenova event. 

Even though eo is exponentially sensitive to N^. this uncertainty does not pose a big problem 
for estimating the gravitational wave signal at its maximum. Indeed, even if we set eo = ei, the 
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Figure 12: The contour plot of constant gravitational wave signal from axion transitions between 
the 6g and the 5g levels for a black hole located at 20 Mpc away from the Earth. The projected 
sensitivity curves assume 10'^ seconds of a coherent integration time. 



transition terms in (66) are still too small to compete with the superradiance upto ei ~ 10~ , 
so that this uncertainty does not affect the dynamics. Also, to estimate the maximal possible 
signal let us concentrate on the very last episode of spinning down which is terminated because 
the black hole reaches the / = 4 Regge trajectory (and not by the Bosenova event, as happens 
for earlier episodes of spindown). At the end of this episode superradiance shuts down and only 



the transition terms in (66) are left. At this point the eo/ei ratio is small, but it starts growing 
as a result of transitions. The transition signal reaches its maximum when eo ~ ei ~ 10~^ and 
then decreases because the occupation fraction ei for the 6g-leve\ drops down. The duration of 
the signal at peak intensity is determined by the transition rate and is of order 

which is of order a day for a stellar mass black hole with a ~ 1 and eo ~ 10~^. 
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In Fig. [12] we present the contour plot for the gravitational wave amplitude as determined 
by (62) for different axion masses and values of a (equivalently, for different black hole masses), 
assuming ei = eo = 10""^ and taking 20 Mpc as a distance to the source. We presented also the sen- 
sitivity curves of various planned gravitational wave detectors assuming the coherent integration 
time 10^ seconds. 

We see that future experiments will be sensitive to the transition signal from the superradiant 
cloud over a large range of axion and black hole masses. What is especially exciting is that the 
Advanced LIGO detector which is scheduled to start operating around 2014 will be probing the 
heavy mass regime for axions, in particular the QCD axion. We will discuss the Advanced LIGO 



reach for the QCD axion in more details in section [44 

Finally, as we already said, as a consequence of Bosenova, every transition between Regge tra- 
jectories goes through a series of 10-100 spin-down episodes interrupted by silent intervals needed 
to build up the cloud. Depending on the distance to the source one may see also gravitational 
wave signal from the earlier episodes, although to study this possibility requires a more detailed 
analysis of the dynamics (in particular, accurate prediction of the e-folding number A^^e). If these 
signals can be observed, then by measuring the frequency and the amplitude of the signal as well 
as the duration of the active and silent intervals one may hope to extract not only the mass of the 
axion, but to estimate the scale as well. 

Also the Bosenova event by itself gives rise to a gravitational wave burst. Assuming that the 
collapse of the cloud happens on a time scale of order Vg the power emitted in gravitational waves 
during the Bosenova event can be estimated as 

Pbn ~ Give'M^^r;^ , 
which translates into the gravitational wave amplitude at the Earth of order 



r VlO-V \WMq) \ r ) 



with the frequency being of order ^ . This signal may be observable for supermassive black holes. 
To conclude, let us emphasize that in estimating the signal strength we were using the free 



particle approximation for axions all the way up to e ~ 10^^, when the bound (51) gets saturated 
and this approximation may be not adequate, as we already discussed. This may be especially 
important for the transition signal where one needs to follow several levels simultaneously. Conse- 
quently, our encouraging estimates here should be considered as a strong motivation for a further 
careful numerical analysis of the system, rather than as an accurate prediction for the signal. 



4.3 Direct observation of the cloud 

Another potential observational consequence of superradiance is the presence of the cloud itself, 
which could be directly detected by precision mapping of the near horizon black hole metric. Such 
a mapping will be made possible by future low frequency gravitational wave detectors, such as 
LISA or AGIS, during the last stages of the inspiral of a compact object (black hole/neutron 
star/white dwarf) into a supermassive black hole. With LISA sensitivities, it will be possible to 
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observe hundreds of such events per year for different galaxies, and to trace up to ~ 10^ orbits in 
an individual event. This will allow mass and spin determination with 10~^ accuracy, and about 
6-7 higher multipole moments of the metric can be measured with better than a few percent 
precision. 

In principle, it is straightforward to calculate the modification to the waveform of the inspiral 
signal due to the presence of the axion cloud. In the regime when non-linearities can be neglected 
(and this is the only regime, where the cloud can stay for a cosmologically long time) the shape of 
the cloud is determined by the well-known wave functions of the hydrogen atom. It is likely that 
the best chances to observe the presence of the cloud are for black holes with moderately small 
values of a//. Indeed, at smaller values of a non-linear effects and processes leading to dissipation 
of axions from the cloud get suppressed allowing for longer lifetimes and a larger cloud mass. On 
the other hand, the size of the cloud grows at small a and the total mass becomes smaller for the 
same value of the spin ("the balerina effect"), making it challenging to see the effect of the cloud 
at too small values of a. A dedicated study is required to find the optimal value of a and to see 
whether the effect is observable. 

Another subtlety with using gravitational wave signal from compact inspirals to detect the 
presence of the cloud is that the non-spherical metric perturbation induced by an infalling compact 
object may be strong enough to cause a mixing between superradiant and non-superradiant levels 
and induce the disappearance of the cloud, similar to the effect of axion self-interactions. 

To summarize, directly probing the structure of the cloud with extreme mass ratio inspirals 
is an interesting possibility awaiting for a dedicated theoretical study to decide on whether it is 
feasible. Another possibility worth exploring are the effects of the cloud on the accretion disk of 
stellar mass black holes. In this case, the cloud could excite resonant modes of the accretion disk, 
the so called quasi-periodic oscillations. 



4.4 The QCD axion and superradiance 

The QCD axion is the best motivated of all axion-like particles and by itself serves as one of the 
major motivations for the whole axiverse framework, so let us summarize here what range of its 
parameter space will be probed by on-going and future black hole observations. 

Unlike for other axions, the QCD axion mass [la and decay constant ja are related to each 
other by (|2]). Furthermore, non-perturbative string corrections to the QCD axion potential take 
the form 

y string ^ A^C""^ COS 0/ fa , 

where the energy scale A is either Planck or string scale (in exceptional cases it might be sup- 
pressed by the SUSY breaking scale, A^ ~ MpiFsusY [21]) and the instanton action S in explicit 
constructions is bounded from above as 

S<'^. (67) 

Ja 

S is close to saturating the above bound if a compactification manifold is neither too anisotropic 



nor strongly warped. It was suggested [12] that the upper bound (67) follows from very general 



41 



properties of quantum gravity. However, to the best of our knowledge, there is no specific proposal 



for the exact numerical coefficient that should appear in a conjectured sharp version of (67). 

To solve the strong CP problem the instanton action should be sufficiently large, S > 200. 
Combined with the above arguments this suggests that the scale fa for the QCD axion is unlikely 
to be significantly higher than few x 10^^ GeV or equivalently, that the QCD axion mass is 
unlikely to be significantly lighter than 10~^° eV. Also, from a bottom-up perspective, the values 
of fa close to the grand unification scale, corresponding to masses rua ~ 3 ■ 10~^° eV, appear to 
be well-motivated. 

4.4.1 Black hole spindown and Advanced LIGO 

These arguments motivate thinking of the consequences of superradiance in an axion mass range 
that is as heavy as possible while still affecting black hole dynamics. Of course, from a purely 
phenomenological approach any limit on the QCD axion parameters in the high fa regime are 
still extremely interesting. The current measurements of black hole masses and spins, presented 
in Fig. [3} already suggest an upper bound on the axion decay constant at the level 

fa<2- 10'' GeV . (68) 

For higher values of the decay constant, i. e. lighter axion mass, the gap in the upper panel of 
Fig. [3] would shift towards heavier black hole masses and would contain rapidly spinning black 
holes inside. Of course, at the moment one should consider this bound as indicative. First, the 
data points in Fig. |3] may be subject to significant systematic uncertainties. For instance, the 
highest spin black hole in Fig. |3] is GRS 1915+105, and the model for the soft X-ray spectrum 
of this object suggests a much smaller value of the spin a/vg ^ 0.56 [6], than the one presented 
here (from [5]). Second, rapidly spinning black holes in Fig. |3]may turn out to be young enough 
to stay in the gap region — for instance, the age of the second fastest spinning black hole in Fig. [3] 
(LMC X-1) is quite short — of order 5 million years |40j . 

All these uncertainties will get rectified with more data coming. Given that a black hole 
produced as a result of the stellar collapse can be as light as ~ 2M0, spin measurements alone 



can potentially improve the bound (68) by a factor of few. Still, the above theoretical arguments 
suggest that it may be not enough to discover the QCD axion. A plausible situation could be 
that the QCD axion is light enough to affect the dynamics of the lightest stellar mass black hole 
through superradiance, but is still too heavy to produce a noticeable gap in the spectrum of 
rapidly rotating black holes. This makes it especially important to study other consequences of 
superradiance that may allow to discover the QCD axion in such a situation. 

One exciting possibility to discover the QCD axion is through observing of the gravitational 



wave signal from superradiating black holes at Advanced LIGO, as discussed in section 4.2[ Es 



timates presented there indicate that Advanced LIGO may see the transition signal for the QCD 
axion, although the annihilation signal will probably have too large frequency to be observed at 
that experiment. As we mentioned, details of the transition signal — such as the strength of the 
signal at the maximum intensity — may even provide an estimate for the decay constant, which 
would be a further confirmation that the signal is related to the QCD axion, rather than to some 



other axion-like particle. In Fig. fflwe zoomed in the high mass region of the plot in Fig. 12 relevant 
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for the QCD axion. This plot shows only the estimated signal for the 6g — )■ 5g transitions. For 
heavy axions the higher I levels are likely to be relevant; our estimates show that the corresponding 
signal is very similar to the one in Fig. |4j Note that we pick 10^ second as a coherent integration 
time for this plot, so that the actual reach can be even better for longer integration times. We 
see that Advanced LIGO has the potential to probe the most interesting mass range for the QCD 
axion, and this is even more true for the Einstein telescope. 

Coming to other probes, unfortunately, it is impossible to use gravitational waves to probe 
the shape of the QCD axion cloud around stellar mass black holes as it could be done for lighter 
axions affecting supermassive black holes. However, it is interesting to study whether the QCD 
axion cloud may resonantly excite certain characteristic perturbations in the accretion disc that 
would signal about the presence of the cloud. 



4.4.2 Photon signals: radio waves, 7- and X-rays 

Finally, there could be QCD axion specific signatures related to the direct coupling of the QCD 



axion to Standard Model fields. First, there is a coupling (30 ) to photons. As we already discussed 
this coupling converts axions from the superradiant cloud into photon with the rate given by (31). 
This conversion rate is too slow to affect the dynamics of superradiance, but it may still provide 
an observable signal on Earth. The photons in question are almost monochromatic radio waves. 



and the production rate (31) translates in the following flux at Earth, 



radio 



W f lO'' GeV y /^x f B y /IkpcV 
m2 V fa J V 10- V V 10^ Gauss; \ r J 



There are several challenges for this signal to be observable, and they all arise because the frequency 
of the signal is equal to the axion mass, and corresponds to radio waves with wavelengths of order 
at least few hundred meters. The first difficulty arrises because these wavelengths cannot be 
observed on the Earth's surface, since the ionosphere is not transparent for radio waves at these 
low frequencies. However, this problem may be solved by using space- (or Moon-) based radio 
telescopes. The major challenge, however, is that the electron density around a black hole should 
be quite small, Ue < 100 cm"^, for the signal to be able to escape from the source. It is very 
hard to find a source satisfying this property, given that one needs substantial magnetic field in 
the vicinity of the cloud for efficient axion-to-photon conversion. This possibility will be studied 
in P]. 

Probably a more promising possibility is related to the coupling ([T]) of the QCD axion to 
the QCD instanton number. To observe the consequences of this coupling would be especially 
interesting given that the interaction ([T]) is a genuine footprint of the QCD axion. An intriguing 
possibility to achieve this is related to the following observation. As follows from the estimate 



(47) the axion field in the cloud may reach values of order 



~ 2 J fa , (69) 

i.e., the ratio (j)/ fa becomes of order one. For the QCD axion this ratio is nothing else as but the 
local value of the CP violating 6'-parameter in the QCD. QCD properties are rather different at 
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large values of the 6'-parameter — for instance, the pion mass is smaller by a factor of order 3 at 
6' = TT as compared to the 6^0 vacuum, where we live (see, e.g., jH] for a recent discussion). 
Given that (multi)pion exchange is one of the dominant forces responsible for the nuclear binding 
it is natural to expect that nuclear binding energies change by order one in the regions with ^ ~ 10 
Consequently, it is natural to speculate that some of the stable nuclei may become unstable 
towards disintegration, and gamma- or beta-decay when they enter in the region of the cloud as it 



approaches the maximum size. Even if (47) somewhat overestimates the maximum 9 in the cloud, 
it appears very probable that 6 becomes of order one at least during the Bosenova events. The 
characteristic timescale for the latter is set by the black hole size and is of order 10~^ seconds. 
Consequently, it's only strong and electromagnetic nuclear instabilities that have enough time to 
be important during the Bosenova event. 

Under the optimistic assumption that an order one fraction of nuclei in the vicinity of the black 
hole horizon decays and produces 7-quanta with MeV energies let us estimate the resulting flux 
of photons at the Earth from the Bosenova event. If the black hole accretes with an Eddington 
limited rate, the total amount of matter within a distance of order from the black hole horizon 



can be estimated as MBufg/rE ~ 10^^ GeV~ 1O~^^M0, where te is the Eddington time (33). 
This may give rise to the emission of order 10^^ photons with MeV energies on a time-scale of 
order 10~^ seconds. This is not very much — for a black hole at 10 pc away one would obtain one 
photon per 10 m^ at the Earth. However, there are several ways the signal can be significantly 
stronger. First, the major limiting factor in the above estimate is the amount of matter in the 
vicinity of the black hole horizon. This amount may be roughly 22 orders of magnitude larger 
immediately after the supernova explosion. Of course, this is a violent event providing lots of 
radiation by itself, and also immediately after the explosion the metric perturbation due to the 
surrounding matter is likely to be strong enough to damp the superradiance. However after the 
environment cleans up a bit this may give rise to a signal many orders of magnitude stronger than 
in the Eddington regime. 



Also the signal may last significantly longer if, as the estimate (69) suggests, nuclei may get 
destabilized not only during the Bosenova collapse, but also when the cloud is still stable. Another 
possibility for the signal to last longer is for the Bosenova to produce long-lived axion clumps (such 
as "pulsons" of that would be able to escape from the near-horizon region. 

Finally, rather than directly detecting photons one may look for spectral X-ray lines of exotic 
elements in the vicinity of the black hole, that could have been formed as a result of the nuclear 
disintegration triggered by the large axion fielcfl 

4.5 Accelerated growth of black holes 

Yet another indirect consequence of superradiance is that axions, if present, would accelerate the 
black hole growth. As a toy illustrative example we pick the Eddington saturated thin disc model 
(see [l6] for a review and references). In this model the black hole mass evolves according to 

dM 1 - eMid) M 



dt euia) Te 



+ Msr , (70) 



•^Note, that this effect was not taken into account in the analysis of 
''We thank Steve Kahn for pointing this out. 
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where, the accreting gas is assumed to have zero metalhcity, te is given by (33), and eAf(a) is 



the radiation efficiency for the accrettion — the fraction of the incoming energy that gets radiated 
away in the process of accretion. It depends only on the dimensionless spin-to-mass ratio 



and in Fig. 13 we have shown the coefficient (1 — eM)/^M appearing in (70) as a function of a. 
Finally, the M^r-term accounts for superradiance, and is given by 



-/i, 



(71) 



where the omitted terms are those related to non-linear effects. 

The time evolution of the black hole spin in the thin disc model is determined by the following 
equation 

da s{a) ^ 



dt eM{,a)TE 



(72) 



where the ratio s/eM as a function of a is also shown in Fig. 13 As before, the a^^-term describes 
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Figure 14: The effect of an axion on the black hole growth history for the simplest thin disc 
Eddington limited accretion (lower curves) and for a more realistic model thin disc model of 
Eddington accretion taking into account result of magnetohydrodynamical simulations (details of 
both models can be found in [17]). 



the effects of superradiance and is equal to 



^5^(a-V-2a)r,iV,. (73) 



Fig. 13 now makes the effect of superradiance evident: the solid line there indicates that the 
mass of a slowly rotating black hole accreting in the Eddington regime grows almost an order of 
magnitude faster than the mass of a rapidly rotating black hole. 

The origin of this effect is easy to understand. For a rapidly rotating Kerr black hole the size of 
the last stable orbit is significantly smaller than for a Schwarzschild black hole of the same mass. 
As a result the accreting plasma radiates a larger fraction of its rest mass before falling into a 
rotating black hole. In the Eddington saturated regime the radiation pressure is the main limiting 
factor for the accretion rate — hence, the accretion proceeds faster for a Schwarzschild black hole, 



as the solid line in Fig. 13 shows 



The dashed line in Fig. 13 indicates that a black hole accreting in the Eddington regime 
rapidly spins up and keeps growing with a high value of spin. Instead, if a black hole is affected 
by superradiance it follows the Regge trajectory, where its spin can be significantly lower and, as 
a consequence, the growth rate is much faster. Given that even a single axion affects a large range 
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of black hole masses, superradiance may significantly affect the growth history of supermassive 
black holes, as illustrated in Fig. [I4j 

It is worth noting that quasars hosting ~ 10^ Mq black holes were observed at as high redshifts 
as z = 6.43, and some authors think that these observations may present a challenge for the 
conventional story of black hole growth [17] . Of course, the actual dynamics describing black hole 
evolution is likely to be significantly more complicted with merger events playing a significant 
role (see, e.g., [38] for a recent overview), and it is premature to decide whether any new physics, 
such as axions, is needed to explain the existing observations. However, it is likely that with 
future X-ray and gravitational wave data as well as with a progress in numerical simulations the 
evolution history of supermassive black holes will be understood much better, and new physics 
might eventually be required, especially if quasars with significantly higher redshifts are to be 
discovered. At any rate, even if a conventional astrophysics is able to explain the data, it is useful 
to keep in mind that axions, if present, are able to significantly affect the story. 

5 Bookkeeping of anthropic axions 

Before concluding, let us elaborate on one particularly interesting feature of axions in the mass 
range relevant for the present paper, fia ^ 10~^^ eV and with a high symmetry breaking scale, 
fa ~ Mgut- Namely, the axion abundance relative to baryons is given bjj^ 

where 7 is an order one coefficient different from unity if an axion is heavy enough, so that the 
effective number of degrees of freedom at the onset of its oscillations is different from the current 
value; da = (f){t = 0)/ fa is an initial axion misalignment angle, and 

PiOa) ^ el (75) 

for small 9a (the shape of P{Oa) for general 9a can be found, e.g., in Fig. 4 of [I]). We see that 
axions with masses significantly heavier than ~ 10^^^ eV would produce a contribution to the 
dark matter density larger than the observed value Qcdm ~ 5^b, unless we happened to have an 
atypically small initial misalignment angle 9a- 

As was realized long ago ^], this does not mean that such axions are necessarily in conflict 
with the observed value of ( = flcdm/^b- Indeed, if inflation lasted sufficiently long (and especially 
if there were a period of eternal infaltion in the past) an initial misalignment angle 9a is a dynamical 
parameter that varies in space on scales much longer than the current size of the Universe, so there 
will always exist regions with sufficiently small value of 9a to be in agreement with the observed 
value of (. Still, one may wonder what is the probability for an observer in such a Universe to 
find himself in a region with as small values of ( as we observe. 

For a single QCD axion this question was addressed a number of times in the past |l9l [50] . 
These treatments differ in some details, however they agree that the observed value of ( does not 

^For the QCD axion this formula gets modified due to the temperature dependence of the axion mass. This is 
not important for the discussion below. 
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appear to be anomalously small. Let us see how this conclusion changes if there are more than one 
axions with masses > 10~^^ eV. Note that depending on which of the parameters are allowed to 
vary the answer may be more or less sensitive to the uknown details of the statistics of string vacua 
and to the infamous ambiguities with the probability measure in an eternally inflating Universe. 

We find the approach of [50] — to keep all the parameters apart from 9a fixed — the safest from 
this point of view. In other words, we are restricting to comparing observers with all microphysical 
parameters the same as ours, but the inflationary dynamics automatically produce different initial 
values of 9a for them. In the string landscape the axion abundance is the last parameter that 
may vary, so this approach is maximally close to the logic applied for predicting the results of a 
conventional lab experiment — we fix all particles physics parameters to the known values and see 
what the dynamics of the system gives us. The important difference with a lab experiment is that 
now we cannot ignore the selection effects — the formation of observers is impossible in the regions 
where the dark matter-to-baryon ratio C, is either too big or too small. These so called anthropic 
boundaries were estimated in [IH]. Namely, for < 2.5 perturbations at the scales close to our 
galaxy's cease to grow, while at C > 100 the density of baryons becomes so small that the disc 
fragmentation instability leading to star formation does not develop. The observed value C ~ 5 
appears to be somewhat too close to the lower end of this interval. 

Unfortunately, to quantify whether there is a real problem, we are still left with an ambiguity 
related to the choice of the inflationary measure — the problem of comparing numbers of observers 
measuring different values of C given that these numbers are infinite in an infinite Universe for 
any C, in the anthropically allowed region. Following [50], let us consider what happens with 
one particular choice — the causal diamond measure of [ST]. This choice amounts to counting 
the number of observers in a single Hubble patch of the late time de Sitter evolution. Another 
simplifying assumption of [50] is that the number of observers per baryon is approximately constant 
for 2.5 < C < 100 and zero otherwise. 

The nice feature of the axion setup is that the prior probability distribution for C, is known. 



Indeed, generalizing (74) to the case when more than one axion is present we find 

C(^a) = 5^c(/x,)P(^,). (76) 

a 

The initial values 9a is getting set during inflation when the axion backreaction on the cosmological 
expansion is negligible, so that the prior distributions for all 9a are flat. Then the probability to 
observe the dark matter-to-baryon ratio smaller than the observed C = 5 value is equal to 

V = U--! -fl^, (77) 

where a factor (1 + C(6'a))~^ is specific to the causal diamond measure and appears because the 
total number of baryons within a horizon at the transition to the de Sitter regime is proportional 
to this factor. The normalization factor M is equal to 



2.5<C(6»a)<100 



1 + aoa) 
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Expressions (77), (78) are significantly simplified in the limit when all anthropic axions are suf- 



ficiently heavy, so that the approximation (75) is accurate. In this regime one can get rid of the 



mass dependence in (76) by rescaling 9a — )■ c{fia) ^^'^6 a- As a result, after integration over angular 



variables in the 9a space, one obtains, 

nn) = 'oo = 0.3, 0.16, 0.06, 0.02, 0.006, . . . (79) 

J2.5 i+c 

where n is the number of axions, and we presented the numerical value of the probability for the 
first few values of n. We see that the probability drops exponentially as the number of axions in 
the anthropic window grows, however, remains high enough for the first few values of n. Clearly, 
this general trend — that at large number of axions the probability distribution is peaked at the 
higher anthropic boundary for ( is generic and independent of the choice of the measure. It is 
just a consequence of a geometrical factor (^("~2)/2 numerators of integrals in (79). For 



instance, if we droped the (1 + C(^a))~^ factor and just used the prior probability distribution for 
C (restricted to the anthropically allowed region), we would get 0.08, 0.03, 0.007, 0.001 for the 
first few probabilities. It is worth pointing out that these probabilities are sensitive to the position 
of the anthropic boundary at large (, which is not the case in the presence of a single axion, as 
pointed out in [50] . 

If a sufficiently large number of axions is to be discovered in the anthropic region or, even if, 
for a single QCD axion, it turns out that a significant fraction of cold dark matter is composed 
of WIMP's, these probabilities may start being problematically low. However, one should keep 
in mind that there are lots of uncertainties in the above estimates. Apart from a choice of the 
infiationary measure, the assumption that the number of observers per baryon is constant over 
the whole anthropic interval appears to be a vast oversimplification, due to both astrophysical and 
astrobiological reasons. On the astrophysical side it is far from clear that the number of stars is 
proportional to the number of baryons in the whole range 2.5 < C < 100- Furthermore, the number 
of observers may not scale linearly with the number of stars both due to astrophysical reasons, for 
example due to close encounters, and due to astrobiological, if the early stages of the evolution 
of life can be significantly accelerated by the possibility of the transfer of organic molecules (or 
primitive forms of life) from one stellar system to another (given that the closest known planetary 
system is just 10 light years away this possibility is neither necessarily hypothetical nor untestable). 

To summarize, we see that at the current stage of affairs there is no reason to be discouraged on 
the possibility of a discovery of multiple anthropic axions with astrophysical black holes observa- 
tions. Conversely, if several anthropic axions were to be discovered (or even a single one if WIMPs 
constitute a significant fraction of dark matter) this will provide us with serious motivation to 
scrutinize how the number density of observers depends on the baryon-to-dark matter ratio. 

6 Conclusions 

We hope to have convinced the reader that black hole superradiance for axions is an extremely 
rich phenomenon that has good chances to be observed in near future measurements of black hole 
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properties. Ongoing black hole spin measurements may trace gaps or Regge trajectories in the 
spectrum of rapidly spinning black holes. Advanced LIGO may observe gravitational wave signals 
from the QCD axion cloud around stellar mass black holes as far as the Virgo cluster for masses 
down to 10"^" eV which correspond to an axion decay constant close to the grand unification 
scale. In a more distant future, gravitational waves may be observed for supermassive black holes 
at lower frequencies by experiments such as LISA and AGIS. The low frequency gravitational 
wave detectors may also see the effect of the cloud on the waveforms during extreme mass ratio 
inspirals. Finally, for the QCD axion the superradiant cloud might also give rise to direct photon 
signals. 

In this paper our main goal has been to develop a general intuition about superradiance 
development and its consequences without going into an extensive numerical work. Given the 
richness of the system it seems inevitable that detailed numerical simulations will be required in the 
future to obtain accurate quantitative predictions. This is especially important for predicting the 
strength, duration and precise waveform of gravitational wave signals from superradiant clouds. 
Simulations are needed both for the accurate prediction of superradiance rates, which to large 
extent has already been accomplished in |30j, and most importantly to get an accurate description 
of the superradiant cloud including axion self-interactions. It's worth stressing, however, that some 
of our result are very robust. For instance, the black hole Regge trajectories of Fig. |3]are mostly 
determined by the basic superradiance condition ([s]) and do not depend on the above uncertainties 
(apart from the left most declining segments of these curves; however, in that region even larger 
uncertainties are likely to come from variations in the accretion rate for different black holes). 

Of course, we expect also other qualitative results obtained here to reproduce well the gross 
features of the system, although given its richness more surprises are possible. Most importantly, 
our results appear encouraging for prospects to observe superradiance with future astrophysi- 
cal data, and this justifies further theoretical efforts for better understanding of this fascinating 
process. 
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Appendix A: Calculating the tunneling exponent (22) 



In principle, it is straightforward to calculate the tunneling integral / numerically, however, let 
us also describe an anlytical method that works for near-extremal black holes, a ^ rg and for the 
frequency u right at the boundary of the superradiance region u = mwj^. 

The latter condition has the following meaning. As we see from Fig. [5} each of the superradiant 
levels has a maximum rate close to the boundary of the superradiant region. Of course, exactly at 



the boundary the width of the level becomes zero. However, it is the pre-factor in (25) that turns 
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Figure 15: The tunneling exponent as a function of the black hole spin. 



zero at u = mw+, while the exponent / just passes smoothly through that point. Consequently, 
by calculating J at a; = mw+ we will find the upper envelope of the family of superradiant rates 
for different levels. Related to this, we will also set I = m, because this corresponds to the fastest 
superradiance rate at any a. Finally, we also set /Iq = u, which given the above assumptions 
corresponds to taking / = 2a. This should be a good approximation given that superradiant levels 
are close to be non-relativistic. Note that for the fastest level the radial number n also grows with 
/, so that the upper bound Q never gets saturated. All these assumptions were also made in [19j . 
The simplification at a; = mw^ is that the location of ri is known, namely ri = r_|_, for this 



choice of parameters. This fact is straightforward to check explicitly using (21). It is also easy to 



understand intuitively — the only way for the tunneling rate (and, consequently, for the imaginary 
part of an eigenfrequency u) to vanish is for the tunneling to be "kinematically forbidden", and 
this is exactly what happens if V{r+) = 0. For all other values of u the potential at the horizon 
is negative, V{r^) = —{u — mwj^Y the imaginary part is non-zero. 



A further simplification happens at a 



-in this case the second turning point coincides 



for Z > 2. This implies that at a; = mw+ one cannot set a 



Tg before 



with the horizon, r2 = 

performing the integral in (22). However, this makes the integration simple. Indeed at w = mw+ 
we can write 

V{r) = (r — r+)(r2 — r)v{r) , (80) 
where in the limit a ^ Vg the function v{r) has a finite non-vanishing limit in the whole interval 
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(r+,r2). Hence at w = mw+, we can write the integral (22) in the hmit a — )■ as 

= 27rt;(r,)i/V2 



I = v{r+Y'\rl + a^) [ ' dr 



r2 — r 1 



r2 — r_ 



r_i_ — r_ 



Now, taking the second derivative of (80) with respect to r and setting a = r = we obtain, 

a{a — 1) 



2r4 

5 



^2) 



Similarly, taking the mixed second derivative of (|80|) with respect to r and a we obtain, 

3a + 1 



ri - a^Oa — + 



1 — a 



(83) 



Equivalently, (83) imphes that at a ~ 1 one has 

3a + 1 



r2 = r„ + 



a — I 



^4) 



Finally, by plugging (82) and (84) into (81), we obtain the following answer for the tunneling 



integral in the extremal Kerr geometry 

I = Tx{2a- v^2a(a - i; 



^5) 



In principle, we can continue as above and work out higher order terms in the (1 — a) -expansion 
to arrive at the approximate analytical WKB formulaes for superradiant rates in the near-extremal 



case. Instead, in Fig. 15 we present the result of a numerical integration of (22) as a function of 
the black hole spin. 
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